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METHODS FOR PREDICTING THE BIOLOGICAL, 
CHEMICAL, AND PHYSICAL PROPERTIES OF 
MOLECULES FROM THEIR SPECTRAL PROPERTIES 

5 

RELATION TO PRIOR APPLICATIONS 

This is a continuation-in-part of application serial no.09/496,314, filed February 
10 1,2000. 

FIELD OF THE INVENTION 

The invention relates to methods for predicting the biological, chemical, and 
physical properties of molecules from their spectral properties. 

15 

BACKGROUND 

The number of known chemical compounds is vast, and the number is 
increasing constantly because methods for isolating and synthesizing molecules 

20 continue to improve. For instance, chemists are now able to employ the techniques of 
combinatorial chemistry to synthesize thousands of different chemical compounds at 
one time using a mixture of only a few interchangeable chemical building blocks. 
Furthermore, chemists are now able to use combinatorial computer models to 
generate large numbers of chemical structures that may in theory be synthesized. 

25 While there are many chemical compounds, only a relative few of those 

compounds may exhibit a particular desirable property, such as pharmaceutical 
activity. Random testing of known compounds to find those that show pharmaceutical 
activity is very expensive and time-consuming. Similarly, there is a need to screen 
compounds for toxicity, so that rational decisions can be made regarding the use and 

30 regulation of compounds that have toxic potential. At present, only a fraction of 
known compounds have been thoroughly tested for their toxicological and potential 
therapeutic properties. To address this problem, scientists have developed methods 
which attempt to predict which compounds are likely to exhibit a particular property. 
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Methods for predicting the properties of chemical compounds are generally 
based upon the related observations that the structure of a compound is related to its 
biological, chemical, and physical properties, and that compounds of similar structure 
exhibit similar properties. These observations have been used to search for new 
5 compounds exhibiting a particular property. For instance, a benzene ring is present 
in both acetaminophen and salicylamide, both of which are analgesics. Although 
incorporating a benzene ring into a new molecule increases the likelihood that it too 
will exhibit analgesic activity, this deduction only narrows the compounds to be 
tested. This approach is still basically one of trial and error, because many 

10 compounds with a benzene ring are not analgesics. Moreover, analgesics without a 
benzene ring will be missed in the search. 

Quantitative structure-property relationships and quantitative structure-activity 
relationships (collectively QSAR) are attempts to quantify the observed relationships 
between the structure of chemical compounds and the extent to which 

15 those compounds exhibit certain properties. For instance, a QSAR might attempt to 
quantify how the analgesic activity of known analgesics that contain a benzene ring 
(such as acetaminophen and salicyamide) depends upon the number and identity of 
substituents on their benzene rings. Once established, such a QSAR could be used to 
predict the analgesic activity of other compounds that contain benzene rings, and 

20 identify those compounds that warrant further investigation as analgesics based on 
their predicted analgesic activity. 

In the terminology associated with the QSAR method, the property for which a 
prediction is sought, such as analgesic activity, is termed the "endpoint." In general 
the endpoint may be any measurable biological, chemical or physical property. 

25 To establish a QSAR, endpoint values are obtained for a set of compounds and 

a correlation is then sought between the endpoint values and some measure(s) of 
structure available for each of the compounds. The measures used to describe or 
reflect the structure of the compounds for which a correlation is sought are termed 
structure descriptors. Structure descriptors may be defined directly with reference to 
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the known structures of the compounds, or may indirectly reflect the structure 
through a property of the molecule that is sensitive to changes in structure. For 
example, an investigator might try to correlate the analgesic activity of compounds 
that contain a benzene ring with either a direct measure of the structure, such as the 
5 number of hydroxyl groups attached to the benzene ring, or an indirect measure of the 
structure of the compounds, like water solubility. If the direct measure is chosen, the 
attempted correlation could only mclude those compounds with hydroxyl groups on 
the benzene ring, while the mdirect measure is more general and could be used to 
include all benzene ring containing compounds in the attempted correlation. 

10 The endpoint data and the structure descriptor (s) for the set of compounds that 

are chosen to establish a QSAR are termed the training set. In general, the reliability 
of a QSAR increases as the number of compounds in the training set increases. If 
possible, the training set desirably includes compounds that exhibit a wide range of 
endpoint values and possess diverse structures. 

15 If the endpoint and structure descriptor(s) are sufficiently correlated, a 

mathematical or graphical representation may be obtained. For example, the growth 
inhibition of certain gram negative bacteria by aromatic amines is correlated in a 
linear fashion with the logarithm of the octanol-water partition coefficient, and the 
correlation indicates that as the aromatic amines become more hydrophobic they are 

20 more likely to inhibit the growth of these bacteria (Hansch and Leo, Exploring QSAR: 
Fundamentals and Applications in Chemistry and Biology, American Chemical 
Society, 1995, p. 416). 

The QSAR representation may be used to predict endpoint values for other 
compounds from their structure descriptor(s), and the reliability of a QSAR may be 

25 tested using a validation set of data. A validation set includes structural and endpoint 
data for compounds that were not part of the training set. The validation set desirably 
exhibits a diversity of endpoint values and structures that is commensurate with the 
training set. The QSAR is tested by how reliably it predicts endpoint data for the 
validation set compounds from the validation set structural data. For example, a- 
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naphthylamine, an aromatic amine, is about two times less toxic than predicted by its 
octanol-water partition coefficient, indicating that the simple linear relationship 
described above is not always reliable and that some specific interaction is responsible 
for its behavior. 

5 The Hammet equation is another example of a simple QSAR that relates a 

single structure descriptor (in this case, derived directly from the structure of the 
compounds) to an endpoint (in this case, the equilibrium or rate constant for a 
particular type of reaction). The Hammet equation is a Imear equation of the form: 

In Kx = pa + In Kh 

10 

The electronic parameter (a) is a measure of the ability of a group of atoms (typically 
a substituent at a particular structural position) to donate (negative value) or withdraw 
(positive value) electron density to or from the reaction center of the molecule to 
which it is attached. The slope parameter (p) is a measure of the sensitivity of the 

15 reaction to the withdrawal or release of electron density, and is constant for a 

particular type of reaction. Kx and Kh are, respectively, the equilibrium constant for 
reaction of the substituted molecule and the equilibrium constant for reaction of the 
unsubstituted parent molecule. The reaction rate constants, kx and kn, may replace 
the equilibrium constants. A plot of In Kx (or In kx) versus a is often linear and may 

20 be used, for example, to predict the equilibrium constant (or rate constant) for other 
structurally similar compounds from only the a values of their substituents. 

The Hammet parameter a is an example of a structure descriptor that is 
derived from experimental data for compounds of known structure. As described 
above, a is a measure of the electronic properties of a substituent. Specifically a 

25 measures the ability of a substituent at a particular position to donate or withdraw 
electron density and may be defined according to the following equation: 

a=logKx-logKH 
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where Kxand Kh are, respectively, the acid ionization constants for an aromatically 
substituted benzoic acid derivative and for unsubstituted benzoic acid. The parameter 
a for amino group substitution in a position para to the reaction center is determined, 
for example, from the measured acid ionization constants of j7-aminobenzoic acid and 
5 benzoic acid using the above equation. Although depending on the position of 
substitution, electron-withdrawing substituents generally tend to stabilize the anion 
formed when benzoic acid ionizes, making Kx larger than Kh and a positive. 
Conversely, electron-donating substituents tend to destabilize the anion and generally 
have negative a values that are also dependent upon the position of substitution. Acid 
10 dissociation data for benzoic acid derivatives with various substituents in ortho, meta 
or para positions have been used to generate a values for many substituents in these 
structural positions. Unfortunately, the a values derived for substituents in this 
manner are typically not valid for multiply substituted molecules because a values 
may not be additive. 

15 Due to the inability of a values to accurately reflect the cumulative electronic 

effect of multiple substituents, alternative structure descriptors have been investigated. 
For example, Ijzerman et al. derived structure descriptors from the assigned ^^C 
NMR shifts of the aromatic carbons in a family of N-r^rr-butylphenylethanolamines 
(Ijzerman et al., 7, Med. Chem. , 29: 549-554, 1986). Ijzerman et al. subtracted the 

20 value of the chemical shift for the carbon atoms in benzene (128.5ppm) from the 
assigned chemical shift values for each of the carbons in the benzene ring of each of 
the substituted N-?^rr-butylphenylethanolamines to yield sets of structure descriptors 
for the compounds. A correlation was found between certain linear combinations of 
these structure descriptors and the p2-andrenoreceptor intrinsic activity of the 

25 substituted N-f^r^-butylphenylethanolammes. Since these structure descriptors are 
defined in a manner similar to the Hammet a parameter, they necessarily rely upon 
knowledge of structure beforehand. Also, the scope of usefulness for this type of 
structure descriptor is Innited to the compounds containing a benzene ring. 
Furthermore, these structure descriptors ignore the aliphatic carbon atoms of the 
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molecules, and their contribution to the intrinsic activity (a valid assumption only 
when comparing the intrinsic activities of compounds that differ minimally in 
structure). Additionally, the correlations of Ijzerman et al. have been questioned 
because the structure descriptors they derived from the chemical shifts of the aromatic 
5 carbons in these compounds are themselves highly correlated, and their use leads to 
the unrealistic result that substitution at one ortho position of the benzene ring affects 
intrinsic activity in a manner opposite to substitution at the other ortho position (see, 
Cramer et aL, Quant. Struct Met. Relat., 7: 18-25, 1988). 

A similar attempt to utilize spectroscopically derived parameters as structure 

10 descriptors was made by Nishikawa and Tori (Nishikawa and Tori, /. Med, Chem. , 
27: 1657-1663, 1984), The reactivity of the P-lactam ring of 3-substituted and 3- 
methylene substituted cephalosporins toward alkaline hydrolysis showed a correlation 
with changes in the ^^C NMR chemical shift of selected carbon atoms, and with 
changes in the infrared (IR) stretching frequency observed for the P-lactam carbon- 

15 oxygen bond. Selection of these measures was guided by knowledge of the structure 
of these compounds, and the role p-lactam reactivity plays in their antibiotic activity. 
Like the method of Ijzerman et al. , the method of Nishikawa defines its structure 
descriptors in terms based upon a reference compound of known structure and utilizes 
spectral features assigned to particular atoms of known structures. Furthermore, the 

20 method, through its definition of the structure descriptors, is similarly restricted to 
narrow classes of compounds that differ only by simple substitutions. 

Structure descriptors may be obtained theoretically or experhnentally. Dipole 
moments and lowest unoccupied molecular orbital (LUMO) energies are examples of 
structure descriptors that may be obtained from theoretical quantum mechanical 

25 calculations on known structures. Experimental data correlated with a specific 

structural feature, common to a set of closely related compounds, may also be used to 
generate measures of structure (for example, the Hammett a parameter). Bulk 
experimental measures of structure, such as partition coefficients (as a measure of 
polarity) and molar refractivities (as a measure of steric size) can also be utilized as 
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Structure descriptors. Structure descriptors based upon bulk physical properties have 
the advantage that they do not require structural knowledge beforehand, hovs^ever such 
descriptors lack specificity. For example, compounds with vastly different biological 
activities may have very similar partition coefficients. 
5 A particularly important type of biological QSAR uses, as an endpoint, the 

ability of one molecule to bind to another molecule. An example of such an endpoint 
would be the ability of a series of molecules to act as ligands for a regulatory protein, 
such as a hormone receptor. In recognition of the important role played by three- 
dimensional structure, especially for biochemical reactions where molecular 

10 recognition is an important factor, the field of 3D-QSAR has emerged. The 3D- 
QSAR technique is exemplified by the Comparative Molecular Field Analysis 
(CoMFA) method of Cramer and Wold (U.S. Patent No. 5,025,388). The CoMFA 
method attempts to correlate the three-dimensional steric and electrostatic properties 
of a series of molecules with their relative endpoint values. The steric and 

15 electrostatic properties of a molecule are obtained from quantum mechanical or 
electrostatic calculations based upon known molecular structures and serve as 
structure descriptors. The calculations, in effect, map the electron density 
distribution around a molecule to create a 3-D picture of its steric and electrostatic 
fields (collectively, the molecular field). Those steric features (e.g., bulky 

20 substituents) and/or electrostatic properties (e.g., a strong molecular dipole) that are 
most important in determining endpoint values are revealed by comparing the 
molecular fields of the molecules in the training set to their endpoints. The advantage 
of calculated 3D-QSAR molecular field structure descriptors is that unlike structure 
descriptors referenced to a certain structural feature, molecular field structure 

25 descriptors enable the identification of structurally dissunilar molecules that have 
similar steric and electrostatic properties. 

However, a particular problem associated with the CoMFA method and other 
3D-QSAR techniques is that these methods generally require some assumptions about 
how molecules orient themselves relative to each other upon binding. Selecting the 
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common alignment of a training set containing diverse structures may be problematic, 
leading to incorrect predictions of binding ability. Furthermore, QSAR based upon 
quantum mechanical or electrostatic potential calculations also suffers to some extent 
from the inaccuracy of the calculations themselves. These calculations are by nature 
5 approximate, and become less and less reliable as molecular size increases. In 

addition, the effects of solvation on the quantum mechanical properties of a molecule 
are often difficult to calculate and typically are ignored. Finally, these calculations 
are time-consuming and require knowledge of the molecular structures beforehand. 
Isolation and structure elucidation of molecules is an expensive and time- 

10 consuming process that is desirably avoided, especially when screening large numbers 
of molecules (e.g., those generated in a combinatorial library). Furthermore, many 
molecular structures are proprietary, so there may be reluctance on the part of 
investigators to share training set structural and endpoint data with competitors in a 
common effort to improve the reliability of predictive QSARs. 

15 It would be advantageous to have QSAR methods that utilize structure 

descriptors which may be obtained without knowledge of or disclosure of molecular 
structures. Such structure descriptors could be more specific than bulk physical 
properties, and enable the QSAR methods to differentiate dissimilar molecular 
structures that exhibit similar bulk physical properties. Furthermore, the structure 

20 descriptors could avoid definition with reference to a particular structural element 
(e.g., the Hammet a parameter), and therefore be usefiil for establishing QSARs for 
more structurally diverse sets of compounds. Additionally, such QSAR methods 
could eliminate the costly step of structure elucidation and obviate the need for 
isolation of the subject compounds. 

25 QSAR methods that utilize structure descriptors that are inherently reflective 

of the steric and electrostatic properties of molecules and the effects of solvation 
thereon are also needed. Such methods would also desirably eliminate the necessity 
for assumptions regarding molecular orientation in relationship to intermolecular 
binding and obviate the need to rely upon approximate calculations. 
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SUMMARY OF THE DISCLOSURE 

The present invention avoids some of the foregoing problems by providing a 
method for predicting a biological activity of a molecule, by obtaining spectral data 
(such as NMR data) for a test compound, and comparing the spectral data for a test 
5 compound to a pattern derived not exclusively from the assigned spectral data (such 
as NMR data) of a training set of compounds having known biological activity. 
Similarities between the pattern of spectral data associated with the biological activity 
of the training set compounds and the spectral data for the test compound are detected 
to determine whether the test compound is predicted to exhibit the biological activity, 

10 The spectral data of the compound for which a prediction is sought need not first be 
correlated with corresponding structural features. Furthermore, the pattern of 
spectral data associated with the biological activity may be derived without first 
correlating the spectral data with corresponding structural features. Training set 
patterns and similarities between the training set patterns and the test compound's 

15 spectral data are conveniently detected, in some embodiments, by segmenting the 
spectral data of the training set and test compounds into sub-spectral units (bins). 
These sub-spectral units, or bins, may be of a width corresponding to the digital 
resolution of the method used to generate the spectral data or greater. 

For example, the biological activity of a test compound may be predicted by 

20 comparing the signals in bins of the training set spectral data that are found to be 
associated with a biological activity (such as strong estrogen receptor binding) to 
signals in corresponding bins of the spectral data of the test compound. Numerous 
signals in the test compound's spectral data that fall within bins corresponding to 
strong estrogen receptor binding are an indication that the test compound possesses 

25 strong estrogen receptor binding. Some of the bins of the training set spectral data 
may contain signals more consistently associated with strong estrogen receptor 
binding, in which case the presence of signals in the corresponding bins for the test 
compound would be more heavily weighted in assigning a predicted biological activity 
to the test compound. Lesser numbers of signals in the test compound's spectral data 
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that fall within bins corresponding to strong estrogen receptor binding indicate that 
the test compound possess only moderate or weak estrogen receptor binding. 

The spectral data of the training compounds and the test compound may be 
just one type of spectral data (such as NMR, for example ^^C-NMR), or more than 
5 one type of spectral data (such as a composite of two or more of NMR, mass spectral, 
infrared, ultraviolet- visible, fluorescence, or phosphorescence data). In particular 
embodunents, the spectral data is a composite of two or more of nuclear magnetic 
resonance spectroscopic (NMR) data, mass spectroscopic (MS) data, infrared (IR) 
spectroscopic data, and ultraviolet-visible (UV-Vis) spectroscopic data. 

10 In other embodiments, the spectral data of the training set compounds is 

segmented into sub-spectral units (bins), and scaled to normalize the importance of 
different signals (e.g., those from different types of spectra in a composite or those 
arising from signals of different inherent intensities) prior to pattern recognition. In 
particular embodiments, the scaling is auto-scaling. In yet other embodiments, the 

15 spectral data of the training set compounds is weighted prior to pattern recognition to 
emphasize those sub- spectral units (bins) that are most important for differentiating 
endpoint classes (such as strong versus weak estrogen receptor binding) of compounds 
in the training set. In more particular embodiments, the weighting is Fisher- 
weighting. 

20 A pattern of spectral data associated with a biological activity can be 

advantageously derived from the training set spectral data using computer 
implemented pattern recognition techniques. Furthermore detection of similarities 
between the pattern derived from the training set spectral data and the spectral data 
exhibited by a test compound is also advantageously performed using computer 

25 implemented methods. In one such approach, the pattern of spectral data associated 
with a biological activity is derived for a training set of compounds by segmenting the 
spectral data and generating a set of canonical variate factors, one for each bin of the 
segmented spectral data. These canonical variate factors are used with the spectral 
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data of a test compound to yield a prediction of the biological activity of the test 
compounds. 

The methods of the present invention are advantageously computer 
implemented. In one such embodiment, the method is a computer implemented 

5 system for predicting iiological activity of a test compound, in which input spectral 
data is received for a test compound, and for a set of training compounds having a 
known biological activity. The spectral pattern derived from the training set (derived, 
for example, using computer implemented pattern recognition programs) and the 
spectral pattern of the test compound are compared to determine whether the spectral 

10 patterns of the test compound match spectral patterns of the training set associated 
with a biological activity. The spectral data for the test compound and the spectral 
data for the training set may conveniently be divided into substantially identical 
spectral bins, so that a signal within individual corresponding spectral bins is 
compared between the pattern derived from the training set and the test compound. 

15 In this computer implemented system, the spectral patterns are obtained by 

inputting spectral data such as one or more of nuclear magnetic resonance data, mass 
spectral data, infrared data, ultraviolet-visible data, fluorescence data, and 
phosphorescence data. In particular embodiments, the spectral data is a composite of 
nuclear magnetic resonance data and mass spectral data. The spectral data of the 

20 training set is converted into principal components (PCs) and canonical variates 
(CVs). Peaks in particular bins of the canonical variates that are associated with a 
biological activity (such as high affinity to a hormone receptor) are identified, and the 
test compound is analyzed for the presence of one or more (and ideally many) 
corresponding peaks in its composite spectrum. 

25 Another aspect of the invention is a method for predicting a biological, 

chemical, or physical property of molecules, by providing spectral data segmented 
into spectral sub-units, for a plurality of training compounds; inputting the segmented 
spectral data and endpoint data into a pattern-recognition program; training the 
pattern-recognition program with the segmented spectral data and endpoint data to 
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establish a relationship between the spectral sub-units of the segmented spectral data 
and the endpoint; providing segmented spectral data for a test compound that is 
segmented into substantially the same spectral sub-units that were used for the 
training set, and comparing the relationship between the spectral sub-units of the 
5 segmented spectral data and the endpoint to the spectral sub-units of the test 
compound's segmented data to predict the endpoint of the test compound. In this 
method, the structures of the training compounds and the test compound are not 
necessarily known beforehand. In particular embodiments of the method, the 
segmented spectral data of the training set is auto-scaled and Fisher-weighted before 

10 training the pattern recognition program. 

The spectral data is chosen from the group consisting of nuclear magnetic 
resonance data, mass spectral data, infrared data, UV-Vis data, fluorescence data, 
phosphorescence data, and composites thereof, for example ^^C NMR data, EI MS 
data, and composites thereof. The endpoint may be a ligand-target molecule-binding 

15 affinity, such as estrogen-receptor binding affinity. Other examples of the endpoint 
are a measure of biodegradability; a measure of toxicity; participation in a metabolic 
pathway; a partition coefficient; a reaction rate; a quantum yield, a measure of 
phototoxicity, an equilibrium constant; and a site of reaction on a molecular structure. 
In a particular embodiment, the endpoint is the octanol/water partition coefficient. 

20 Additional non-spectral structure descriptors may be utilized along with 

segmented spectral data to provide an expanded set of structure descriptors useful for 
establishing a predictive relationship for the endpoint. Examples of non-spectral 
structure descriptors that do not necessarily require structural knowledge beforehand 
include partition coefficients, solubilities, relative acidities, relative basicities, pKa, 

25 pKb, reaction rates, and equilibrium constants. In a particular embodiment, the 

partition coefficient is the octanol/water partition coefficient. Calculated non-spectral 
descriptors are one example of descriptors requiring structural knowledge beforehand 
that may be utilized along with segmented spectral data in establishing a predictive 
relationship for an endpoint. 
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Another aspect of the invention is a method of using spectral data as a set of 
structure descriptors for a compound that does not necessarily require knowledge of 
the compound's structure beforehand, by providing spectral data of a training set of 
compounds and segmenting the spectral data into bins. 
5 Yet another aspect of the invention is a method for establishing a relationship 

between spectral data and a biological, chemical, or physical property, by providing 
spectral data for a training set of compounds, segmenting the spectral data into bins, 
and detecting patterns in the bins of the spectral data that are associated with the 
property. The method may also include detecting correspondmg patterns in spectral 

10 data of test compounds to select the test compounds having the property. 

Furthermore, the test compounds may be mixtures of compounds. The spectral data 
of the training set can be auto-scaled and weighted (for example by Fisher-weighting) 
to help identify data that are most strongly associated with the biological activity. 
Knowledge of the structural features that lead to the spectral data is not needed 

15 beforehand. 

Yet another aspect of the invention is a method of determining the structural 
features of a plurality of compounds that contribute to determining a particular 
endpoint property exhibited by the compounds, by providing segmented spectral data 
for the plurality of compounds; providing endpoint data for the plurality of 

20 compounds; establishing a spectral data-activity relationship (SDAR) by identifying 
the segmented spectral features that bias toward increased endpoint values and the 
segmented spectral features that bias toward decreased endpoint values; and 
identifying the structural feature leading to the segmented spectral features that bias 
toward increased or decreased endpoint values for the plurality of compounds. As 

25 with the earlier described embodiments of the invention, the segmented spectral data 
may be a composite of several types of spectral data. 

Any of the methods of the present invention can be performed without 
reference to a chemical structure of the test compound. Hence the spectral features of 
the training set and the test compound may be compared, without determining the 
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chemical structure of the compounds of either the training set or the test compound or 
compounds. Any of the foregoing methods can also be incorporated into a computer 
readable medium, having stored thereon instructions for performing the steps of these 
methods. 

5 The foregoing and other objects, features, and advantages of the invention will 

become more apparent from the foUowuig drawings and the detailed description of 
several embodiments. 

BRIEF DESCRIPTION OF THE DRAWINGS 

10 Figs 1(a), 1(b) and 1(c) are the ^^C NMR, EI MS, and IR spectra, repectively, 

for bisphenol A (4,4'-isopropylidenediphenol, structure shown in Fig. 1(a)). 

Figs. 2(a), 2(b) and 2(c) show tables of data that correspond respectively to 
each of the spectra shown in Figs. 1(a), 1(b) and 1(c). 

Fig. 3 shows a hypothetical set of structure descriptors derived from the 
15 spectral data summarized in the tables of Figs. 2(a), 2(b) and 2(c). 

Fig. 4 shows a flowchart of a particular embodiment of the spectral data- 
activity relationship (SDAR) method. 

Fig, 5 shows the discriminant function using ^^C Nuclear Magnetic Resonance 
(NMR) spectral data for 30 compounds in the Spectral Data- Activity Relationship 
20 (SDAR) model. The X-axis is the first canonical variate (discriminant function) and 
the Y-axis is the component frequency. The numbers in each box correspond to the 
numbers in Table 1 that serve to identify the compounds used in the SDAR. White 
boxes correspond to strong estrogen receptor binding compounds, and gray boxes 
correspond to medium estrogen receptor binding compounds. 
25 Fig. 6 presents the first canonical variate factor weights using ^^C NMR 

spectral data for 30 compounds in the SDAR model. The X-axis is the bin number 
and the Y-axis is the factor weight relative intensity. The bins are numbered from 
550 (corresponding to 0 ppm) to 770 (corresponding to 220 ppm). 
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Fig. 7 shows the discriminant function using composite ^^C NMR spectral data 
and Electron Impact (EI) mass spectral data for 30 compounds in the SDAR model. 
The X-axis is the first canonical variate and the Y-axis is the component frequency. 
The numbers in each box correspond to the numbers in Table 1 that serve to identify 
5 the compounds used in the SDAR. White boxes correspond to strong estrogen 

receptor binding compounds, and gray boxes correspond to medium estrogen receptor 
binding compounds. 

Figs. 8(a) - 8(b) present the first canonical variate factor weights using 
composite ^^C NMR spectral data and EI mass spectral data for 30 compounds in the 

10 SDAR model. The X-axis is the bin number and the Y-axis is the factor weight 

relative intensity. EI mass spectral data are in the bins numbered from m/z 0 to 550. 
^^C NMR spectral data occupies the bins numbered from 550 (corresponding to 0 
ppm) to 770 (corresponding to 220 ppm). 

Fig. 9 shows the discriminant function using ^^C Nuclear Magnetic Resonance 

15 (NMR) spectral data for 107 compounds in the Spectral Data-Activity Relationship 
(SDAR) model. The X-axis is the first canonical variate (discriminant fiinction) and 
the Y-axis is the second canonical variate. The symbol S represents a strong estrogen 
receptor binder, M represents a medium estrogen receptor binder, and W represents a 
weak estrogen receptor binder. 

20 Fig. 10 presents the first canonical variate factor weights using ^^C NMR 

spectral data for 107 compounds in the SDAR model. The X-axis is the bin number 
and the Y-axis is the factor weight relative intensity. The bins are numbered from 
550 (corresponding to 0 ppm) to 770 (corresponding to 220 ppm). 

Fig. 11 shows the discriminant function using composite ^^C NMR spectral 

25 data and Electron Impact (EI) mass spectral data for 107 compounds in the SDAR 
model. The X-axis is the first canonical variate and the Y-axis is the second 
canonical variate. The symbol S represents a strong estrogen receptor binder, M 
represents a medium estrogen receptor binder, and W represents a weak estrogen 
receptor binder. 

15 
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Figs. 12(a) - 12(b) present the first canonical variate factor weights using 
composite ^^C NMR spectral data and EI mass spectral data for 107 compounds in the 
SDAR model. The X-axis is the bin number and the Y-axis is the factor weight 
relative intensity. EI mass spectral data are in the bins numbered from m/z 0 to 550. 
5 ^^C NMR spectral data occupies the bins numbered from 550 (corresponding to 0 
ppm) to 770 (corresponding to 220 ppm). 

Fig. 13 is a diagram of a distributed computing environment in which the 
present invention can be implemented. 

Fig. 14 is a block diagram of a computer system that can be used to implement 
10 the present invention. 



DETAILED DESCRIPTION OF 
SEVERAL DISCLOSED EMBODIMENTS 



15 The illustrative embodiments that follow are best understood by first 

understanding terms in the art that are relevant. These terms are defined below. 



Definitions and Abbreviations 

20 LOO - leave-one-out cross validation 

RBA - relative binding affinity 

NMR - nuclear magnetic resonance 

IR- mfrared (spectroscopy) 

MS-mass spectrometry 
25 EI MS -electron impact mass spectrometry 

UV-Vis -ultraviolet-visible (spectroscopy) 

CODESSA - comprehensive descriptors for structural and statistical analysis 
SDAR -spectral data-activity relationship 
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Endpoint - a particular biological, chemical, or physical property or a set of 
such properties for a compound that are either qualitatively or quantitatively 
measurable. 

Structure Descriptors - any direct or indirect measure of the structure of a 
5 compound that may be obtained by theoretical or experimental means. 

Training Set - endpoint data and structure descriptors (structural data) for a 
group of compounds used to establish a correlation between the endpoint property and 
the structures of the compounds. 

Validation Set - endpoint data and structure descriptors (structural data) for a 
10 group of compounds used to test the reliability of the correlation between the endpoint 
property and the structures of the compounds. 

Segmented Spectral Data - spectral data that is divided into discrete sub- 
spectral units (bins), each of which spans a particular spectral range. The spectral 
range spanned by a particular bin corresponds to a range of frequencies or a range of 
15 wavelengths for spectroscopic data and may be equal to the digital resolution of the 
spectral data or greater. For mass spectrometric techniques, the spectral range within 
each bin corresponds to a particular mass or range of masses and may be equal to the 
digital resolution of the spectral data or greater. The bins need not all be of equal 
width. The spectral data that is divided into bins may either encompass all the 
20 spectral data of a particular type that are available or cover only a portion of the 
spectral data of a particular type that are available. Each bin contains information 
derived from the spectral signals (or lack thereof) that appear within the spectral 
range defined by a particular bin. The structural aspect(s) of the compounds that give 
rise to the information falling within any particular bin need not be known. 
25 Spectral Data-Activity Relationship (SDAR) - a correlation between 

endpoint data and spectral data for a group of compounds, useful for among other 
things for predicting the endpoint data for compounds from their spectral data. 

Nuclear Magnetic Resonance (NMR) - a phenomenon exhibited by a large 
number of atomic nuclei in which nuclei in a magnetic field absorb energy from a 
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radio-frequency field at certain characteristic frequencies. Particular examples of 
nuclei that exhibit this phenomenon include ^^C, and ^^P. 

Mass Spectrometry (MS) - a method of chemical analysis in which the 
substance to be analyzed is placed in a vacuum and reduced to low pressure. The 
5 resulting vapor is exposed, for example, to a beam of electrons which causes 
ionization to occur, either of the molecules or their fragments. The ions thus 
produced are accelerated and then passed through a mass analyzer that separates the 
ions according to their mass. 

Electron Impact Mass Spectrometry (EI MS) - a mass spectrometric 
10 technique in which the ionization of molecules and their fragments is accomplished by 
a beam of electrons that mipacts the molecules and their fragments. Typically, as the 
energy of the electron beam is increased, the number of fragments produced from a 
molecule increases. 

Infrared Spectroscopy (IR) - an analytical technique which measures a range 
15 of wavelengths (or frequencies) in the infrared region or near-infrared region of the 
electromagnetic spectrum that are absorbed by a specimen, which characterize its 
molecular constitution. Infrared absorption bands identify molecular structure 
components, such as aromatic, olefin, aliphatic, aldehyde, ketone, carboxylic acid, 
alcohol, amine, and amide groups. The frequency at which absorption occurs also 
20 reflects the frequency at which the bonds m these components stretch and and/or 
bend. 

Ultraviolet-Visible Spectroscopy (UV-Vis) - an analytical technique which 
measures a range of wavelengths (of frequencies) in the ultraviolet and visible regions 
of the electromagnetic spectrum that are absorbed by a specimen, which characterize 
25 the electronic energy levels of its molecular constituents. UV-Vis absorption bands 
may be characteristic of certain molecular components, such as aromatic groups or 
carboxyl (CO) groups. 

Fluorescence Spectroscopy - an analytical technique which measures a range 
of wavelengths (or frequencies) of light a molecule emits in passing from a higher to 
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lower energy electronic state during a given time period (such as the first millisecond) 
after absorbing a photon of light. Fluorescence wavelengths and emission intensity 
reflect the redistribution of energy in the molecule after light absorption. 
Fluorescence excitation spectroscopy reflects the efficiency with which a molecule 
5 converts absorbed energy into fluorescent emission as a function of the wavelength of 
the absorbed photons. 

Phosphorescence Spectroscopy - an analytical technique which measures a 
range of wavelengths (or frequencies) of light a molecule emits in passing from a 
higher to lower energy electronic state on a time scale beyond the first millisecond 

10 after absorbing a photon of light. Phosphorescence wavelengths and emission 
intensity also reflect the redistribution of energy in the molecule after light 
absorption. Phosphorescence excitation spectra reflect the efficiency with which a 
molecule converts absorbed energy into phosphorescent emission as a function of the 
wavelength of the absorbed photons. 

15 Principal Component Analysis (PC A) - this pattern recognition technique 

and the related technique of partial least squares analysis reduce the structure 
descriptors used in SDAR to a smaller number of linear combinations of the original 
set of descriptors that are called principal components (PCs). The optimum number 
of principal components will yield the smallest standard error of prediction while 

20 capturing as much of the variance shown by the structural data as possible. (See 
generally, Kramer, R., Chemometric Techniques for Quantitative Analysis, Marcel 
Dekker, Inc., 1998) 

Comprehensive Descriptors for Structural and Statistical Analysis 
(CODESSA) - a set of various structural descriptors typically utilized in three 

25 dimensional quantitative structure-activity studies that includes constitutional 

descriptors, topological descriptors, geometrical descriptors, electrostatic descriptors, 
and quantum mechanical descriptors, all of which require knowledge of structure 
beforehand. (See , Tong et aL, /. Med. Chem., 39: 380-387, 1995 and CoUantes et 
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al., y. AnaL Chem, , 68: 2038-2043, 1996, both of which are incorporated herein by 
reference) 

Auto-scaling - a method whereby the quantitative spectral information 
contained within each particular bin is compared for all compounds in the training set 
5 to yield an average value and a standard deviation. Then, for each bin comprising the 
structure descriptors of a given compound, the quantitative spectral information 
therein is expressed as a number of standard deviations above or below the average 
for each bin. Autoscaling equalizes the importance of inherently weak spectral 
signals falling within certain bins with the importance of inherently strong spectral 

10 signals falling within certain other bins in describing a set of spectrally derived 
structure descriptors. It may also equalize the importance of different types of 
spectral data in a composite of spectral data. 

Fisher-weighting - a method whereby the quantitative spectral information in 
bins that are important for classifying the training set compounds into different 

15 endpoint groups, such as strong and medium binders to the estrogen receptor, are 

enhanced. For each bin, the variance of the quantitative spectral information between 
the endpoint groups is divided by the variance of the quantitative spectral information 
within the endpoint groups. The resulting dividend becomes a weighting factor that 
has a magnitude larger than one when a particular bin has an important role in 

20 distinguishing the endpoint groups. Each bin is multiplied by its weighting factor to 
yield structure descriptors that are more sensitive to subtle but significant spectral 
variations. 

Leave-one-out (LOO) Cross- Validation - a method whereby each compound 
in the training set is systematically excluded from the data set, after which its 
25 endpoint value is predicted by the spectral data-activity relationship derived from the 
remaining compounds (See, Cramer et al.. Quant, Struct-Act, Relat, 7: 18-25, 1998, 
incorporated herein by reference). Cross-validation is useful for judging the 
reliability of a spectral data-activity relationship, especially where a validation set of 
compounds is not available. 

20 
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EXAMPLES 

The following examples further illustrate the QSAR methods of the present 
invention. The methods utilize spectral data as structure descriptors and correlate the 
5 spectral data with specific biological, chemical, or physical endpoints without the 
need to assign spectral features to their corresponding structural elements. In 
recognition of the advantages presented by using spectral data as structure descriptors, 
a correlation provided by the methods described herein is termed a Spectral Data 
Activity Relationship (SDAR). 

10 In various embodiments, experimental and/or calculated molecular spectral 

data may be used as structure descriptors. If experimental spectral data is used in a 
particular embodiment, there is no need to actually know the molecular structures 
beforehand. Calculated molecular spectral data may be utilized as a surrogate for 
experimental data in some embodiments. In other embodiments, calculated molecular 

15 spectral data may be useful, in combination with an SDAR derived from experimental 
spectral data, for screening compounds that are created as part of a computer 
generated combinatorial library. In yet other embodiments, calculated molecular 
spectral data may be generated for compounds for which no experimental spectral 
data is available and used to establish an SDAR or to screen those compounds for a 

20 particular property. 

In some embodiments, spectral data is used as a set of descriptors, for 
example descriptors of molecular structure. The pattern of the spectrum is 
determined, for example by segmenting the spectral data into portions covering 
particular spectral regions (e.g., ranges of frequency, wavelength, chemical shift, 

25 mass to charge ratio, and the like). The number and/or the intensity of the spectral 
signals within each segmented region may serve as the structure descriptors. The 
spectral data may be of one type or be a composite of several types, including NMR, 
MS, IR, UV-Vis, Fluorescence, and Phosphorescence. Spectral data of a particular 
type may be utilized in its entirety or in part. While spectral data are often used to 



21 



WDN:AAR:th 7/31/2000 4239-55412 NIH (DHHS) Ref. No. E-209-99/0 Express MaU No. EL#121363839US 

Date of Deposit: July 31, 2000 

elucidate the structure of the compound that yields them, the information contained in 
the spectra may be used in some embodiments without the need to interpret the 
spectra. Furthermore, the spectral data may be used in certain embodiments without 
the need to know the structures of the compounds beforehand. Segmented spectral 
5 data is particularly amenable to encryption for secure analysis. 

In some embodiments, ^^C NMR, EI MS, and IR spectra may be combined to 
yield a set of structure descriptors. Examples of these types of spectra for the 
compound bisphenol A are shown in Figs. 1(a), 1(b), and 1(c) respectively. Each 
type of spectrum contributes information regarding the structure and electronic 

10 properties of a molecule. For example, in the ^^C NMR spectrum shown in Fig. 1(a), 
the peaks reflect the number of different types of carbon atoms and the electronic 
environment of each type of carbon atom in bisphenol A. Peaks exhibiting the largest 
5 (chemical shift) values correspond to carbon atoms of a type that are the least 
shielded from the applied magnetic field by the electrons of the molecule. The EI MS 

15 spectrum shown in Fig. 1(b) reveals the molecular mass (m/z = 228, 228 amu) of 
bisphenol A and the mass of the fragments into which the molecule breaks up when 
the electrons impacting the molecule deliver enough energy to break bonds. For 
example, the most prominent peak in the mass spectrum shown in Fig. 1 (b) appears 
at a m/z of 213. The peak at m/z=213 corresponds to the fragment formed by the 

20 loss of a methyl group (CH3, m/z = 15) from bisphenol A (i.e. 228-213 = 15). The 
mass spectral peaks thus reflect the fragments of structure that comprise the molecule. 

The IR spectrum shown in Fig. 1(c) reveals the various types of bonds that are 
present in bisphenol A. For example the peak that appears at 3368 cm"^ in the IR 
spectrum corresponds to the two oxygen-hydrogen bonds found in bisphenol A. 

25 While the ^^C NMR, EI MS and IR spectral features have been discussed with 
reference to the structure of bisphenol A to illustrate the types of structural 
information contained in these spectra, it should be noted that in some embodiments 
the structure of the molecule responsible for the features need not be known and no 
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attempt need be made to assign particular spectral features to corresponding structural 
features. 

The spectral data from Figs. 1(a), 1(b) and 1(c) are compiled in tabular form 
in Figs. 2(a), 2(b) and 2(c), respectively, as tables of spectral frequencies and relative 
5 intensities. Such compilations are typical of the output of modern digital 

spectrometers. The data as presented in these tables may be used to construct a set of 
molecular structure descriptors for bisphenol A and similar sets of data may be used 
to generate structure descriptors for any number of compounds. 

To construct the set of structure descriptors shown in Fig. 3, each type of 

10 spectrum is segmented into spectral "bins" and the information content of each bin 
becomes a separate structure descriptor. Fig. 3 shows a set of raw structure 
descriptors for bisphenol A derived, for illustrative purposes only, from the spectral 
data as presented in the tables of Figs 2(a), 2(b) and 2(c). It is possible to mclude 
more information from the spectra than considered here, such as the intensity of IR 

15 absorption at frequencies that do not correspond to a distinct peak. Each bin in this 
case corresponds to a particular integer m/z ratio for the EI MS data, a range of 
chemical shift frequencies for the ^^C NMR data, and a range of wavenumber 
frequencies for the IR data. Moreover, each bin contains the relative intensity of the 
spectral feature(s) fallmg within a range covered by a particular bin. In this example, 

20 bms 1 to 49 were not used (because such low mass ions are common and may also 
result from air contamination), bins 50 to 549 correspond to m/z ratios from 50 to 
549, bins 550 to 770 correspond to 1 ppm segments of the ^^C NMR spectrum over 
the range of 0 to 221 ppm, and bins 771 to 1120 correspond to 10 cm"^ segments of 
the IR spectrum over the range from 4000 cm'^ to 500 cm"\ Those spectral bins that 

25 do not contain spectral data are given a value of zero and are not shown in Fig. 3. 
Other regions of the spectra can be selected to correspond to the bins in other 
embodiments of the invention, which is not limited to the specific examples recited 
herein. In addition, the spectra can be divided into different segments, instead of the 
integer segments that are merely given as examples. 
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Tracking specific spectral data from the various spectra to the set of structure 
descriptors in Fig, 3 helps to illustrate how, in some embodiments, the spectral data is 
utilized to construct the structure descriptors. For example, the intensity of the peak 
at m/z=213 in the mass spectrum appears in bin number 213. The intensity of the 
5 peak appearing at 154.2 ppm in the ^^C NMR spectrum is in bin number 704 because 
the peak falls within the range of 154 ppm to 155 ppm which corresponds to bm 
number 704. The peak appearing at 3368 cm'^ in the IR spectrum falls within the 
frequency range 3370 cm^ to 3360 cm^ and thus the intensity of the peak appears in 
bin number 834 since this bin corresponds to this particular frequency range. These 

10 three entries are highlighted and shown in bold in Fig. 3, 

The segmented spectral data structure descriptor set for bisphenol A shown in 
Fig. 3 may be altered in some embodiments to normalize the data from each type of 
spectrum to yield structure descriptors of similar magnitude from each type of 
spectral data. For example, in order to normalize the IR data and ^^C NMR data to 

15 the greatest value in the MS data (100), they might be expressed as the number of 
peaks falling within a particular bin's range, multiplied by 100. The 21 appearing in 
bin 834 would thus be replaced with the value 100 as would the 411 appearing in bin 
704, because in both cases only one peak appears in the spectral ranges corresponding 
to these bins. 

20 In other embodiments, UV-Vis, fluorescence, and phosphorescence spectral 

data is segmented into bins alone or in combination with each other and with NMR, 
MS, andIR data. 

In addition to not requiring structures beforehand, experimental spectral data 
as used in some embodiments offers several advantages over calculated spectral data. 
25 Most experimental spectral data reflects the quantum mechanical properties of 

molecules as determined by both structure and solvation. However, mass spectra do 
not reflect solvation. Experimental spectral data also avoids the difficulty and 
approximations inherent in deriving structure descriptors from quantum mechanical or 
electrostatic potential calculations, which are especially inaccurate when the 

24 
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calculations attempt to include solvation effects. Experimental spectral data may be 
reflective of the solution conformations that are responsible for a particular endpoint 
property as well as the role molecular flexibility plays in determining the endpoint 
property. Moreover, experimental spectral data reflect specific structural features and 
5 are therefore preferable to bulk structure descriptors such as partition coefficients that 
lack structural specificity. Additionally, experimental spectral data is already 
available for many compounds in various databases. 

In some embodiments, segmented spectral data is further pre-treated before 
being combined with endpoint data and subjected to statistical pattern recognition or 

10 artificial intelligence based pattern recognition to extract an SDAR. Fig. 4 

summarizes, in flowchart form, a particular embodiment of the steps that may be 
taken to establish an SDAR from spectral data. These steps and others may be 
performed using a computer system. 

In Fig. 4, spectral data is obtained at 20, and then the spectral data is 

15 segmented into bins at 22, which in some embodiments is carried out in a manner 
similar to that described above with reference to Figures 1,2, and 3. The process of 
segmentation is repeated for the spectral data of each of the compounds in the training 
set. The spectral data for each of the compounds of the training set can be segmented 
into a matrix of bins that represent the training set. These segmented spectral data 

20 sets 24 comprise the structure descriptors for the training set compounds. Patterns in 
the matrix of bins can then be determined, for example by pattern recognition 
software. 

In some embodiments, the segmented spectral data 24 is pretreated by auto- 
scaling 26 prior to pattern recognition, to produce auto-scaled, segmented spectral 
25 data 28. Auto-scaling is a method whereby the quantitative spectral information 
contained within each particular bin is averaged over all the compounds in the 
training set to yield an average value and a standard deviation for each bin. Then, for 
each bin comprising the structure descriptors of a given compound, the quantitative 
spectral information therein is expressed as the number of standard deviations above 
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or below the average value for that bin. This reduces the variation in numerical 
magnitude between the bins in the segmented spectral data. Auto-scaling thus helps 
equalize the importance of inherently weak spectral signals that fall within the spectral 
range of certain bins with the importance of inherently strong spectral signals that fall 
5 within the spectral range of other bins. 

In some embodiments, the segmented spectral data is further pre-treated by 
Fisher-weighting 30 the auto-scaled data 28 to improve the ability of pattern 
recognition algorithms to discern which bins (and the spectral information 
corresponding to the spectral range covered by the bin) are most important for 

10 classifying compounds into two or more endpoint classes (such as strong and weak 
estrogen receptor binders). An endpoint class is assigned to each compound based 
upon its endpoint value. For example, compounds could be classified into two 
endpoint classes, those compounds that have an endpoint value above a certain 
number and those that have an endpoint value equal to or below a certain number. 

15 For each bin of the segmented spectral data 24 the variance in value between the 
endpoint classes for the bin is divided by the variance in value within the endpoint 
classes for the bin to yield a weighting factor. The weighting factor has a magnitude 
larger than one when a particular bin has an important role in distinguishing the 
endpoint classes. 

20 Fisher-weighting 30 quantifies the tendency for certain bins (and the spectral 

range corresponding to the bin) to be more helpful than others in deciding the 
endpoint class of a compound. For example, consider a particular bin in which the 
value in the bin is always large for the compounds belonging to one endpoint class 
(such as strong estrogen receptor binders), and always small for the compounds of a 

25 second endpoint class (such as weak estrogen receptor binders). Knowing that the 
spectral data for a compound yields a large value in that particular bin helps to 
classify the compound into the first endpoint class. If on the other hand, the value 
found in the bin is always large for compounds of both classes, it reveals nothing 
about the endpoint class of the compound to know that it has a large value in the bin. 
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In the first case, the variance between the classes would be larger than the variance 

within each class and the weighting factor will be greater than one. Conversely, in 

the second case, the variance between the classes will be about equal to the variance 

within the class, making the weighting factor equal to about one. Each bin of the 
5 segmented spectral data is then multiplied by its weighting factor to yield a set of 

structure descriptors 32 that emphasize the bins most important for deciding the class 

of a compound. Such data is more easily treated by pattern-recognition analysis. 
Pattern recognition 34 is used to establish an SDAR 36, by correlating the 

segmented (and optionally pretreated) spectral data with the endpoint data for the 
10 compounds in the training set. In particular embodiments, the compounds are 
1 1 classified into two or more endpoint classes (e.g. strong versus weak estrogen 

receptor binders) according to their relative endpoint values. The pattern recognition 
y : then determines any segmented spectral features (i.e., the bins) that are characteristic 

of the compounds falling into each of the classes. If Fisher-weightmg is performed 
C 15 prior to pattern recognition, the endpoint classification scheme used for Fisher- 

yj weighting may be retained during pattern recognition. A simplified example is that a 

bin with a strong signal in all the compounds of the test set that have strong estrogen 
Q receptor binding would be a bin that would be considered as predictive of strong 

estrogen receptor binding if a test compound also exhibited a strong signal in that bin. 
20 Multiple such signals would, however, usually be taken into account by the pattern 

recognition software when determining spectral patterns that are associated with the 

endpoint. 

Pattern recognition to establish an SDAR 36 may be accomplished by 
statistical methods or artificial intelligence methods. Underlying these methods is the 
25 idea that if a particular bin corresponds to a spectral signal exhibited by only the 
compounds of a particular class (e.g. strong estrogen receptor binders), the bin will 
bias strongly toward that class. The extent to which a particular bin biases toward a 
particular class depends upon whether or not compounds from all classes (e.g. strong 
and weak estrogen receptor binders) exhibit spectral signals corresponding to the bin. 
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For example, if half of the compounds in one class show a signal that falls within a 
particular bin and slightly more than half of the compounds in another class show a 
signal in the same bin, the bin may only bias slightly toward the latter class. 
Statistical and artificial intelligence methods attempt to quantify these biases and 
5 provide a basis for classifying compounds to an endpoint class according to their 
segmented spectral data. 

Once quantified, the relationship between the segmented spectral data and the 
endpoint may be visualized in various ways, depending upon the particular software 
package utilized for pattern recognition. One way to visualize the extent to which 

10 individual bins of segmented spectral data bias toward endpoint classes is with a 
canonical variate factor plot. An example of a canonical variate factor plot is shown 
in Fig. 6, which is discussed fiirther with respect to Example 1. The SDAR model 
used to generate Fig, 6 was based upon two endpoint classes. The peaks pointing in 
one direction bias toward one class (strong estrogen receptor binding) and the peaks 

15 pointing in the other direction bias toward the other class (moderate estrogen receptor 
binding). The length of the peaks in Fig. 6 correspond to how strongly the peaks bias 
toward a particular endpoint class. For example, the peak appearing in bin 585 of 
Fig. 6 strongly biases toward the class of strong binders while the peak appearing in 
bin 579 of Fig, 6 biases strongly toward the other class of moderate binders. 

20 An SDAR, once established by pattern recognition, may be used in some 

embodiments to quantitatively predict the endpoint value and endpoint class of any 
compound according to its segmented spectral data. In a particular embodiment, 
segmented spectral data of the type utilized to establish the SDAR is considered along 
with the canonical variate function to yield a prediction of endpoint class. For 

25 example, in reference to Fig. 6, a compound may be qualitatively predicted to belong 
to the class indicated by the upward pointing peaks if it exhibits a large number of 
spectral features falling into the upward biasing bins. The SDAR quantifies the 
prediction by considering not only the number of spectral features that fall into 



28 



WDN:AAR:th 7/31/2000 4239-55412 NIH (DHHS) Ref. No. E-209-99/0 Express Mail No. EL#121363839US 

Date of Deposit: July 31, 2000 

Upward or downward biasing bins, but the strength with which each bin biases toward 
a particular endpoint class. 

One example of a way to visualize the ability of an SDAR to correctly classify 
compounds into their endpoint class is with a discriminant function plot. An example 
5 of a discriminant function plot is given in Figure 5, which is discussed in more detail 
with respect to Example 1 . This discriminant function plot is also based upon the 
same two class SDAR used to generate the canonical variate function of Fig. 6. The 
grey and white boxes correspond to the two endpoint classes (e.g. strong and 
moderate estrogen receptor binding) into which the training set data was divided. The 

10 aggregation of grey boxes and the aggregation of the white boxes, as well as the 
separation of the two aggregates, illustrate the ability of the SDAR to discriminate 
between the two endpoint classes. If the SDAR had not found a correlation between 
the segmented spectral data and the endpoint classes, the discriminant function would 
have exhibited a greater mixing of the grey and white boxes. 

15 Unlike methods using spectral data that is referenced to a parent compound's 

spectral data, the methods of the present invention can in certain embodiments reveal 
the importance of particular structural features in determining the endpoint value 
without requiring knowledge of the identity of those structural features beforehand. 
For example, if a particular bin is identified by the canonical variates as showing a 

20 bias toward strong anti-tumor activity, subsequent elucidation of the structure 
responsible for the signal occupying that particular bin is possible. Furthermore, 
spectral data as used in some embodiments is well suited for examining structure 
activity relationships for compounds of diverse structures where several different 
groupings of atoms may elicit similar quantum mechanical features that are important 

25 to a particular endpoint property. 

Another benefit of using spectral data to establish an SDAR according to the 
methods of the present invention is the flexibility it offers. Once spectral data are 
gathered for a set of compounds, an SDAR may be generated with reference to a 
multitude of biological, chemical, or physical endpoints for which data is available. 
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For instance, the same set of spectrally derived structure descriptors for a set of 
compounds may be utilized with toxicity data and antibacterial activity data for these 
compounds to establish two separate SDARs for the compounds, one for toxicity and 
the other for antibacterial activity, 
5 Because the experimentally based SDAR method does not require input of a 

compound's structure, the method yields accurate results without computer modeling 
of a compound's steric or electrostatic properties. The method also removes the need 
to break the compound into secondary structural motifs or assume a particular 
alignment for how a set of molecules will bind to another molecule. 

10 Methods according to the invention additionally provide rapid and inexpensive 

approaches to screen compounds for selected activities or properties. Screening using 
these methods does not require that time and effort be expended in trying to elucidate 
the structure of compounds beforehand, like screening using 3D-QSAR. Screening 
according to the methods of the present invention, by measuring a compound's 

15 spectral data and inputting it into an SDAR, is much more rapid than screening with 
bioassays. Furthermore, once obtained, the same spectral data for a compound may 
be utilized in multiple SDARs to predict various properties of the molecule. 
Additionally, the methods may also be used to screen mixtures or fractions of 
compounds for spectral features associated with a particular property, thereby 

20 obviating the need to spend time and money isolating compounds that show no 

promise of having the particular property. The methods may therefore be extremely 
useful, for example, in combination with experhnentally generated combinatorial 
libraries of compounds. 

When implemented in a real-time web-based computer environment, the 

25 methods of the present invention may provide investigators with a means for rapidly 
and securely estimatmg the molecular activities of proprietary molecular structures. 
Spectral data as utilized in the present invention may be generated for any number of 
compounds and submitted by an investigator, either directly or as a segmented set of 
descriptors, to a central location containing the SDAR software. The spectral data or 
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segmented spectral data may be encrypted by the user and securely transmitted to the 
central location. Because of the speed at which pattern recognition software can make 
predictions based upon an SDAR, a submission of spectral data and retrieval of a 
prediction based on an SDAR can take place in real time over the worldwide web 
5 without the spectral data being stored at the central location. Alternatively, the 

SDAR and the pattern recognition software may be implemented on a single computer 
in an investigator's laboratory. In either case, there is no need to reveal the structure 
of a proprietary compound to obtain a prediction. The ability to provide spectral data 
in a pre-segmented and encrypted form that is difficult to use for structure elucidation 

10 purposes (i.e., industrial espionage) may also encourage investigators to submit their 
spectral and endpoint data permanently to a central location for use as part of larger 
training sets that may lead to more reliable SDARs. 

It will be understood, however, that many variations and modifications of the 
methods can be made while remaining within the scope and spirit of the present 

15 invention. For example, segmentation of spectral data into bins need not occur, but is 
only one example of how spectral patterns may be observed and compared. Other 
approaches for comparing spectral patterns could also be used, for example artificial 
neural networks that may not require gross segmentation of the spectral data per se 
(image analysis). 

20 

Example 1- Use of ^^C NMR and EI Mass Spectrometric Data to 
Produce a Predictive Model of Estrogen Receptor Binding Activity. 

The Food Quality Protection Act, passed in 1996, mandates that the 

Environmental Protection Agency (EPA) develop screening and testing procedures for 

25 endocrine disrupting chemicals (EDCs). An EDC is defined as "an exogenous agent 

that interferes with the production, release, transport, metabolism, binding, action, or 

elimination of natural hormones in the body responsible for the maintenance of 

homeostasis and the regulation of developmental processes." Estrogenic compounds 

represent a significant subset of the EDCs to be tested. Many of these compounds 
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can be screened by determining how strongly the compounds bind to estrogen 
receptors. Thus there is a need to develop inexpensive and rapid methods to screen 
thousands of potential compounds. SDAR according to the invention meets this need. 
The ^^C nuclear magnetic resonance (NMR) spectrum of a compound contains 
5 frequencies that correspond directly to the quantum mechanical properties of the 
molecule and depend largely on the electrostatic features and geometry, including the 
stereochemical configuration, of the molecule. Furthermore, a solution spectrum, 
like NMR, inherently reflects the effects of solvation on the quantum mechanical 
properties. Electron-impact mass spectral (EI MS) data provide a mass-size 

10 description of molecular substructures (and possibly the whole molecule) as well as 
information about the strength of bonds between the atoms of the molecules. The ^^C 
NMR and EI mass spectral data thus represent sets of quantum mechanical structure 
descriptors that reflect the electrostatic and steric properties of a molecule. Such 
experimental data are readily attainable and often already available. 

15 Two Spectral Data-Activity Relationship (SDAR) models were established, 

based on ^^C NMR and EI MS spectra, for 30 estrogenic chemicals whose relative 
binding affinities (RBA) for the alpha (ERa) and beta (ERp) estrogen receptors are 
available. The SDAR models classified the 30 compounds into two endpoint classes, 
those with strong binding affinities for the estrogen receptors and those with medium 

20 binding affinities for the estrogen receptors. The SDAR based on ^^C NMR data 
alone yielded a leave-one-out cross-validation of 90% and the SDAR based on a 
composite of ^^C NMR and EI MS data yielded a cross-validation of 83 % . Two 
compounds that were classified incorrectly in both SDAR models were in a transition 
zone between these classifications. Both SDAR models correctly classified 25 

25 compounds. The models provided rapid, reliable, and simple ways to predict whether 
a compound binds to the estrogen receptor. 
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Procedure 

The estrogenic relative binding affinities (RBAs) of the 30 compounds were 
obtained from previous publications (Kuiper et al., Endocrinology 138:863-870, 
1997; and Tong et al., Endocrinology, 138:4022-4025, 1997). Most of the '^C NMR 
5 spectrometric and EI mass spectrometric data were obtained from the Integrated 
Spectral Data Base System for Organic Compounds web site 
www.aist.go.jp/RIODB/SDBS, the Aldrich Library of^^C and FT NMR Spectra 
(Pouchert and Behnke, Eds., Aldrich Chemical Company, Volumes 1-3, 1993), 
Spectral Data of Steroids (Frenkel and Marsh, eds., Thermodynamics Research 
10 Center: College Station, 1994), and the NIST MS database software version 1.6. 
Experimental ^^C NMR and EI MS data for five compounds were obtained using 
standard methods. 

Specifically, the ^^C nuclear magnetic resonance (NMR) spectral analyses of 
4-hydroxy-estradiol, ICI 164,384, moxestrol, and norethynodrel were performed at 

15 75.46 MHz on a Varian Gemini 300 MHz NMR (Varian Associates, Inc., Palo Alto, 
CA) spectrometer operating at 301 K. The subject compounds were dissolved in 
CDCb or DMSO (dimethyl sulfoxide). The chemical shifts were defined by denoting 
the CDCb peak at 77.0 ppm and the DMSO peak at 39.5 ppm. The spectral width 
was 21,008 Hz with a 2.6-second delay time between acquisitions. The acquisition 

20 time was 0.495 seconds and the number of points acquired was 20,800. 

The samples were also analyzed by direct exposure probe (DEP) mass 
spectrometry (MS). The mass spectrometers were operated hi the electron-unpact 
(EI) mode, with 70 V electron energy. The ion source temperature was set at 150**C. 
Samples, in solution, were applied to the rhenium wire of the DEP and the solvent 

25 was allowed to evaporate before the analysis was begun. MS data were collected 
until the current of the DEP exceeded 500 mA. Norethynodrel, ICI 164,384, 4- 
hydroxy-tomoxifen, and 4-hydroxy-estradiol were analyzed on a Finnigan TSQ 700 
mass spectrometer and moxestrol was analyzed on a Finnigan 4500 mass spectrometer 
(Finnigan Corp,, San Jose CA). 
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The ^^C NMR and El-mass spectral data points were used as comprehensive 
descriptors for structural and statistical analyses (CODESSA) as in Tong et al., /. 
Med. Chem., 39: 380-387, 1995 and CoUantes et aL, 7. Anal Chem,, 68: 2038-2043, 
1996, both of which are incorporated herein by reference). Mass spectrometric data 
5 for mass-to-charge (m/z) ratios from 50 to 549 were used directly as ordered pairs of 
m/z and relative intensity (maximum of 100). Unassigned ID ^^C NMR chemical 
shifts were segmented into bins over a 0 to 221 ppm range. The ^^C NMR 
frequencies were shifted to bins 550 to 770, so bin 550 was the ^^C NMR spectrum 
for frequencies inside 0 to 1 ppm and bin 770 was the ^^C NMR spectrum for 

10 frequencies inside 220 to 221 ppm. 

To save space, the ^^C NMR spectra were saved as sets of ordered pairs each 
consisting of the respective chemical shift frequency in ppm and the respective area 
under the peak. The area under the peak of a specific chemical shift frequency was 
first normalized to an integer. A non-degenerate frequency was assigned an area of 

15 1; a doubly degenerate frequency had an area of 2; and so forth. This was done to 
provide all the spectra with a similar signal-to-noise ratio and to eliminate line-width 
variations due to differences in NMR-instrument field strengths, shimming, 
temperature, pH, and solvents. The bin defined the number of significant and distinct 
chemical-shift peaks inside a ppm range. The optimal range in this example was 

20 found to be 1 ppm to each bin. To normalize the peak intensity of the ^^C NMR 
spectroscopic data relative to the mass spectrometric data, the area of one ^^C NMR 
peak was set to be 25. The number 25 was selected with the objective of scaling the 
maximum value for the ^^C NMR data inside a bin to near 100, which is the 
maximum value for EI MS data. The number of bins used to input the ^^C NMR 

25 spectra was studied by varying the width of the bins from 0.5 to 5.0 ppm. Again, the 
optimal ppm range for this study was found to be 1.0 ppm. Values were normalized 
to a maximum of 100 prior to pattern-recognition analysis. 

The relative binding affinity (RBA) to the estrogen receptor was defined as the 
ratio of the molar concentration of 17-p-estradiol to the competing compound 
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required to decrease the receptor-bound 17-|3-estradiol by 50%, multiplied by 100. 
Thus, 17-p-estradiol has an RBA of 100 and a log (RBA) of 2.0. Strong binders to 
the estrogen receptor were classified as having a log (RBA) over 0.0, and medium 
binders to the estrogen receptor were classified as having a log (RBA) of less than or 
5 equal to 0.0. There were 17 strong binders and 13 medium binders in the training 
set. No weak estrogen-receptor-bmding compounds (i.e., those with a log (RBA) of 
less than or equal to -3) were included. SDAR models that include more compounds, 
including weak-binding compounds, in the training set are better able to accurately 
predict strong estrogen receptor binding (see example 2). 

10 The analysis of the SDAR model was performed by the leave-one-out (LOO) 

cross-validation procedure in which each compound was systematically excluded from 
the training set and its relative binding activity class predicted by the model (Cramer 
et aL, Quant. Struct, -Act. Relat,, 7:18-25, 1988, incorporated by reference). 

The pattern-recognition software used was RESolve Version 1.2 (Colorado 

15 School of Mines, Boulder, CO). The ^^C NMR and EI MS spectroscopic data for all 
30 compounds were input as text files into a computer programmed with the software. 
The spectroscopic data was then auto-scaled and Fisher-weighted prior to principal 
component analysis (PC A). The discriminant analysis was based on canonical variate 
vectors. Leave-one-out (LOO) cross-validation was used to maximize the size of the 

20 training set. 

Autoscaling compared the quantitative response at each mass spectral m/z bin 
or NMR chemical shift bin to all the other bins in the comparison set. An average 
value (with standard deviation) was calculated for each bin. Then, for each spectrum, 
the quantitative response at each bin was expressed as the number of standard 
25 deviations above or below the respective average. This data-pretreatment step 
equalized the weight of consistent variance of signals with inherently small 
magnitudes (25 units for the NMR bin 558 representing a single methyl carbon) to 
those signals with large magnitudes (130,000 area counts at m/z bin 91, possibly 
arising from a tropylium ion fragment). Autoscaling automatically compensates for 
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gross magnitude variations. It was particularly helpful in this example in which two 
completely different types of analytical spectra were formed into a composite set of 
structure descriptors that are representative of molecular characteristics. 

Fisher weighting was used to emphasize those spectral characteristics 
5 important in distinguishing defined endpoint groups. For each mass spectral m/z or 
NMR bin, the variance between groups was divided by the variance within groups. 
The resulting dividend became a weighting factor with a magnitude larger than one 
when the particular m/z or bin has an important role in distinguishing groups. Fisher 
weighting all of the raw spectra before pattern recognition increased the power of 
10 discriminant analysis to classify spectra correctly. In this example, Fisher weighting 
de-emphasized the relative importance of irrelevant spectral information (such as the 
NMR signals from carbon atoms not important to binding to the estrogen receptor). 

Results 

15 Based only on ^^C NMR spectroscopic data, the statistical pattern-recognition 

program with 8 principal components (PCs) had 81. 1 % of the total variance and a 
cross validation of 90%. The mathematical mode of significance was 91,9% at 1, 
with 28 degrees of freedom. The Wilks discruninant criterion was 99.99% . Table 1 
lists the compound names, an identifying number for each compound that is used in 

20 Figs. 5 and 7, the log (RBA) alpha and log (RBA) beta values, the SDAR training 
input endpoint class for each compound ("class mput")(either S for strong binding or 
M for medium binding), the SDAR prediction ("class predicted NMR'') from the ^^C 
NMR data alone and the SDAR prediction ("class predicted NMR and MS") from the 
^^C NMR and EI MS data. Fig. 5 shows the discriminant function for SDAR derived 

25 from the ^^C NMR data alone which illustrates the ability of the SDAR to discriminate 
between the two endpoint classes in the training set. Compounds denoted with a 
white background exhibit strong RBAs, and compounds with a gray background 
exhibit medium RBAs. Fig. 5 shows that the SDAR yielded a large separation 
between the 15 strong (white rectangles on right side) and the 9 medium binders (gray 
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rectangles on the left side); 4 strong and 2 medium RBAs are in a transition zone 
(middle) between strong and medium RBAs. Table 1 also shows that 27 of the 30 
compounds are correctly group-predicted using only, ^^C NMR data. 3|3- 
androstanediol, 2-hydroxy-estradiol and hexestrol are incorrectly predicted to have a 
medium RBA using 8 principal components. Only 3p-androstanediol and hexestrol 
are incorrectly predicted to have a medium RBA using the canonical variate function. 
Strong binding 3p-androstanediol is most likely incorrectly predicted because the 
compound is similar to the medium RBA 3a-androstanedioL Repeating the analysis 
using a larger training set (Example 2) eliminated even this confusion. 
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TABLE 1 



COMPOUND 


Number 


Log(RBA) 


Log(RBA) 


Class 


Class 


Class 






alpha 


beta 


Input 


Predicted 
NMR 


Predicted 
NMR + MS 


Coumestrol 


1 


1.97 


2.27 


S 


S 


S 


3a-Androstanediol 


2 


-1.15 


-0.52 


M 


M 


S 


Sp-Androstanediol 


3 


0.48 


0.85 


S 


M 


M 


4-OH-Estradiol 


4 


1.11 


0.85 


S 


S 


S 


1 7a-Estradiol 


5 


1.76 


1.04 


s 


S 


S 


1 7|3-Estradiol 


6 


2.00 


2.00 


s 


S 


S 


Bisphenol A 


7 


-1.30 


-0.48 


M 


M 


M 


p-Zearanol 


8 


-1.20 


-1.15 


M 


M 


M 


Clomifene 


9 


1.40 


1.08 


S 


S 


S 


2-OH-Estradiol 


10 


0.85 


1.04 


s 


M 


s 


Dehydroepiandroster 


11 


-1.40 


-1.15 


M 


M 


M 


-one 














Diethystibesterol 


12 


2.67 


2.47 


S 


S 


S 


Dienestrol 


13 


2.35 


2.61 


s 


S 


S 


Dihydrotestosterone 


14 


-1.30 


-0.77 


M 


M 


M 


Estriol 


15 


1.15 


1.32 


S 


S 


S 


Estrone 


16 


1.78 


1.57 


S 


S 


S 


Genistein 


17 


0.70 


1.56 


S 


s 


S 


Hexestrol 


18 


2.48 


2.37 


S 


M 


M 


ICI 164,384 


19 


1.93 


2.22 


s 


S 


S 


Methoxychlor 


20 


-2.00 


-0.89 


M 


M 


M 


Moxestrol 


21 


1.63 


0.70 


S 


S 


S 


Nafoxidine 


22 


1.64 


1.20 


S 


S 


S 


Nandrolone 


23 


-2.00 


-0.64 


M 


M 


M 


Norethindrone 


24 


-1.15 


-2.00 


M 


M 


M 


Norethynodrel 


25 


-0.16 


-0.66 


M 


M 


M 


Progesterone 


26 


-3.50 


-3.50 


M 


M 


M 


Tamoxifen 


27 


0.85 


0.78 


S 


S 


S 


Testosterone 


28 


-2.10 


-2.10 


M 


M 


M 


5a-Androstanedione 


29 


<-2.0 


<-2.0 


M 


M 


S 


53-Androstanedione 


30 


<-2.0 


<-2.0 


M 


M 


S 



5 

Fig. 6 shows the canonical variate function for the pattern recognition obtained 
with the ^^C NMR data only. The positive (upward-pointing) peaks in Fig. 6 
correspond to bins that bias toward a strong RBA for binding to the estrogen receptor 
and negative (downwardly-pointing) peaks correspond to bins that bias toward a 
10 medium RBA. In Fig. 6, m/z bins 550 to 770 refer to the ^^C NMR data from 0 ppm 
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to 1 ppm for bin 550 to data from 221 ppm to 222 ppm for bin 770. The canonical 
variate function reveals that the aliphatic CH2 bins 580 to 585 (30 to 35 ppm) have a 
bias toward medium RBA. It also reveals that the methyl CH3 bins 558 and 566 (8 
and 16 ppm, respectively) have a bias toward strong RBA. Many of the benzyl bins 
5 665 to 700 (115 to 150 ppm) have a bias towards strong RBA. Thus these methyl 
groups and benzyl atoms are indicated to be important for strong RBA. 

Based on the composite ^^C NMR and EI MS data, the statistical pattern- 
recognition program with 14 prmcipal components (PCs) included 89.8% of the total 
variance and had a cross validation of 83.3%. The mathematical mode of 

10 significance was 97.6% at 1, with 28 degrees of freedom. The Wilks discriminant 
criterion was 99.89% . Fig. 7 shows the discriminant function for the composite ^^C 
NMR and EI MS data. Compounds with a white background are strong RBAs and 
compounds with a gray background are medium RBAs. Fig. 7 shows a large 
separation between 15 strong (left) and 9 medium (right) binders. There are 2 strong 

15 and 4 medium binders in the transition zone between the strong and medium binders. 
Table 1 shows that 25 of the 30 compounds are correctly predicted using the 
composite ^^C NMR and EI/MS data. Only 3p-androstanediol and hexestrol are 
incorrectly predicted using the canonical variate function. Again, the strong binding 
3p-androstanediol is most likely predicted incorrectly because it is similar to the 

20 medium binding 3a-androstanediol. 

Fig. 8(a) shows the canonical variate fimction used in the pattern recognition 
of the EI MS data. Fig. 8(b) shows the canonical variate function used in the pattern 
recognition of the ^^C NMR data. The negative peaks in Figs 8(a) and 8(b) 
correspond to bins that bias toward a strong relative binder and positive peaks 

25 correspond to bins that bias toward a medium relative binder, to the estrogen 

receptor. In Figs. 8(a) and 8(b), the label "m/z bins 50 to 550" refers to the EI MS 
data, whereas "m/z bins 550 to 770" refers to the ^^C NMR data. In Fig. 8(a), the 
mass canonical variate is split evenly into bms that bias strong and medium binding. 
Many of the canonical variate bins that showed bias in Fig. 2 are also present in Fig. 
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8(b), but they are pointing in the opposite direction. The opposite directions found in 
Fig. 6 and 8(b) are insignificant technically, arising from the pattern-recognition 
program's arbitrary choice of left and right in the corresponding canonical variate 
score plot. 

5 Only 3p-androstanediol and hexestrol are incorrectly predicted in both SDAR 

models. 3p-androstanediol is closest to 3a-androstanediol, 5a-androstanedione, and 
5p-androstanedione, all three of which have medium RBAs. Note however, 3p- 
androstanediol has been inconsistenly identified in the literature as a medium 
estrogen-receptor binder and a strong estrogen-receptor binder (see Miksicek, /. 

10 Steroid, Biochem, Mol. Biol. 49:153-160, 1994 and Kuiper et al., Endocrinology 
138:863-870, 1997, respectively). 

This example illustrates that, using the present method, compounds can be 
classified by biological activity based on similarities in the spectral patterns of the test 
compounds to the spectral patterns of the training set. In particular, additional test 

15 compounds can be classified according to their expected biological activity by 

obtaining the corresponding spectral pattern(s) of the test compound (e.g. NMR and 
MS patterns). Expected biological activity of the test compound can be predicted by 
detecting similarities between the spectral pattern of the test compound and the 
spectral patterns of compounds in the training set that are associated with a known 

20 biological activity. Similarities in patterns of the canonical variates associated with an 
endpoint (e.g. a biological activity such as estrogen receptor binding) derived for the 
training compounds with patterns in the test compound's segmented spectral data, are 
then used to predict whether the test compound will share the endpoint. 

25 Example 2 - Expanded Training Set for Estrogen Receptor Binding SDAR 

An expanded training set of 107 compounds of varying estrogen receptor 
binding affinities was utilized to create an SDAR model from ^^C NMR spectral data 
and to create an SDAR model from a composite of ^^C NMR spectral data and EI MS 
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spectral data. The training set included weak estrogen receptor binders in addition to 
strong and medium estrogen receptor binders. Data for the training set was obtained 
from the following references: Blair et al., Toxicological Sciences, 54: 138-153, 
2000; Hopert et aL, Environmental Health Perspectives, 106: 581-586, 1998; Zava 
5 and Duwe, Nutr. Cane, , 27: 31-40, 1997; and Kuiper et aL, Endocrinology, 138: 
863-870, 1997. 

For this model, strong estrogen receptor binders were defined as having 
log(RBA) > -0.25 and weak estrogen receptor binders were defined as having 
log(RBA) <-3.0. Medium estrogen receptor bmders were defined as having 

10 log(RBA) in between. The compound names, the input class ("Class"), the predicted 
endpoint class based upon ^^C NMR only ("Class predicted NMR"), and the predicted 
endpoint class based upon ^^C NMR plus EI MS ("Class predicted NMR/MS") for the 
training set are listed below in Table 2. Except for the increased number of 
compounds and a change in the mass spectral range utilized, the analysis was 

15 performed in the same manner as the analysis in Example 1. 



Table 2 



Compound 


RBA 


Class 


Class 
predicted 
NMR 


Class 
predicted 
NMR/MS 


1 ,4-Diphenyl-1 ,3-butadiene 




W 


W 


W 


1 ,6-Dimethylnaphthalene 




W 


W 


W 


1 ,8-Octanediol 




W 


M 


M 


2,2'-Dihydroxy4- 
methoxybenzophene 




W 


W 


W 


2,2'-Methylenebis (4- 
chlorophenol) 


-2.45 


M 


M 


M 


2,4,5-T 




W 


W 


W 


2,4-D 




W 


W 


W 


2,4-Dihydroxyben2ophenone 


-2.61 


W 


W 


M 


2-Chloro-4-nnethylphenol 


-3.66 


W 


W 


W 


2-Chlorophenol 


-3.67 


W 


W 


W 


2-Furaldehyde 




W 


W 


W 
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2-OH-Estradiol (3) 




s 


S 


S 


z-rnenyipnenoi 




VV 


\A/ 
VV 


\A/ 
VV 


2-sec-Butylphenol 
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All compounds without a number) Blair, R.M.; et aL ToxJcoIogicaf Sciences, 54:138-153, 2000. 

1) Hopert A-C; et al. Environmentaf Health Perspectives, 106:581-586 1998. 

2) Zava, D.T.; and Duwe, G. /Vafr. Cane, 27:31-40, 1997. 
5 31 Kuiper et a!., Endrocrinology 1 38:863-870, 1 997. 

Based only on ^^C NMR spectroscopic data, the statistical pattern recognition 
program with 27 principal components (PCs) had 93.0% of the total variance and a 
cross validation of 78.5 %. Figure 9 shows the discriminant function for ^^C NMR 
data of 107 compounds. Compounds that are represented by an S are strong estrogen 

10 receptor binders, compounds that are represented by an M are medium estrogen 
receptor binders, and compoxmds that are represented by a W are weak estrogen 
receptor binders. Figure 9 shows a clustering between the 20 strong, 14 medium, and 
73 weak binders. There are cluster overlaps between the weak and strong binders 
and between the weak and medium cluster region. There are a total of five false 

15 positives and three false negatives when using LOO to predict from the 107 
compound SOAR model of estrogen receptor binding usuig only ^^C NMR data. 
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Figure 10 shows the factors associated with the first canonical variate function 
for the pattern recognition of the ^^C NMR data. The positive peaks in Figure 10 
correspond to bins that bias toward a strong RBA for binding to the estrogen receptor 
and negative peaks correspond to bins that bias toward a medium RBA. In Figure 10, 
5 bins 550 to 770 refer to the ^^C NMR data from 0 ppm for bin 550 to 220 ppm for bin 
770. The aliphatic CH2 bins 580 to 585 (30 to 35 ppm) have a bias toward medium 
RBA. The methyl CH3 bins, such as 558 and 566 (8 and 16 ppm, respectively) have a 
bias toward strong RBA. Many of the aromatic bins 665 to 700 (115 to 150 ppm) 
have a bias towards strong RBA. 

10 Based on the composite ^^C NMR and EI/MS data over the range m/z = 100 to 

m/z=549, the statistical pattern recognition program with 27 principal components 
(PCs) included 85.9% of the total variance and had a cross validation of 77.6 %. 
Figure 11 shows the discriminant function for ^^C NMR data for 107 compounds. 
The compounds that are shown by an S exhibit strong RBAs, compounds that are 

15 shown by an M are exhibit medium RBAs, and compounds that are shown by a W 
exhibit weak RBAs. Figure 11 shows a similar clustering between the 20 strong, 14 
medium, and 73 weak binders as seen Figure 9. There are cluster overlaps between 
the weak and strong binders and between the weak and medium cluster region. There 
are a total of seven false positives and three false negatives when using LOO to 

20 predict a compound's endpoint class from the 107 compound SDAR model of 
estrogen receptor binding using ^^C NMR and EI/MS data. 

Figures 12(a) and 12(b) shows the factors associated with the first canonical 
variate function for the pattern recognition of the composite ^^C NMR and EI/MS 
data. The positive peaks in Figures 12(a) and 12(b) correspond to bins that bias 

25 toward a strong RBA for binding to the estrogen receptor and negative peaks 

correspond to bins that bias toward a medium RBA. In figure 12(a), bins 100 to 549 
are the EI/MS data bms and in Figure 12(b), bins 550 to 770 refer to the '^C NMR 
data from 0 ppm for bin 550 to 220 ppm for bin 770. The aliphatic CH2 bins 580 to 
585 (30 to 35 ppm) have a bias toward medium RBA. The methyl CH3 bins, such as 
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558 and 566 (8 and 16 ppm, respectively) have a bias toward strong RBA. Many of 
the aromatic bins 665 to 700 (1 15 to 150 ppm) have a bias towards strong RBA, 

In both 107 compound SDAR models of estrogen receptor binding, a possible 
cutoff for testing experimentally whether a molecule did bind to the estrogen receptor 

5 would be anything with a positive component 1 as depicted in figures 9 and 1 1 . This 
would leave no positives untested in each model. 

A particularly striking result of the analysis using strong, medium, and weak 
estrogen receptor binders in the training set is the good separation between strong 
estrogen receptor binders and the other compounds as seen in Figures 9 and 1 1 . There 

10 is a marked intermixing of the medium and weak estrogen receptor binders in Figures 9 
and 1 1 with respect to the first canonical variate function (X-axis). These results are 
understandable in light of the physical basis for estrogen receptor binding. Estrogen 
receptor binding, like most ligand-receptor binding interactions, is governed by specific 
intermolecular interactions between the ligand and the receptor. The specific 

15 intermolecular interactions are what distinguish strong binders from the other 

compounds. Conversely, there are few or no specific interactions with the receptor on 
which to distinguish the medium and weak estrogen receptor binders from each other. 



Example 3 - Use of ^^C NMR, EI MS, and Infrared Spectral Data to 
20 Produce a Predictive Model of Biodegradability for Monocyclic Chlorobenzene 

Derivatives 



Two recognized factors in the biodegradation of xenobiotic compounds in the 
environment are the physicochemical properties of the compounds themselves, and 
25 the enzymatic activities expressed by the microbial communities in the environment. 
The vast majority of microbial enzymatic reactions may be classified into one of four 
categories: 1) group transfer; 2) oxidation and reduction; 3) elimination, 
isomerization, and rearrangement; and 4) carbon-carbon bond cleavage. The extent 
to which xenobiotic chemicals can serve as substrates for these enzymatic reactions is 
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determined to a large extent by their structure and particularly by the presence of 
bonds that are similar to those found in the natural substrates for the enzymes. 

Infrared (IR) spectroscopic data may be included in an SDAR to provide a set 
of spectrally derived structure descriptors that include the types of bonds present in 
5 the subject molecules. Infrared data also reflect the modes and frequencies of 

vibration that are available to subject molecules. As chemical reactions are often tied 
to particular vibrations of particular bonds, inclusion of the infrared data may 
improve the ability of an SDAR to predict reactivity. 

The biodegradation data for many monocyclic chlorobenzene derivatives in 

10 sediment may be found in the Database for Environmental Fate of Chemicals 

(www.aist.go.jp/RIODB/dbefc). Additional data on the biodegradability is published 
in Biodegradation and Bioaccumulation Data of Existing Chemicals Based on the 
CSCL, Japan Chemical Industry Ecology-Toxicology & Information Center (JETOC), 
Tokyo, Japan, 1992. The half-life period is used as the endpoint for the 

15 establishment of the SDAR and compounds are classified into two endpoint classes as 
readily biodegradable (R) (half-life < 30 days) and not readily biodegradable (NR) 
(half-life >30 days). The endpoint data for 34 chlorobenzene derivatives are given in 
Table 3 below. 

20 Table 3 



Compound 


Half-life (days) 


Input Class 


Monochlorobenzene 


46.2 


NR 


1 ,2-dichlorobenzene 


36.9 


NR 


1 ,3-dichlorobenzene 


433 


NR 


1 ,4-dichlorobenzene 


385 


NR 


1 ,2,3-trichlorobenzene 


23 


R 


1 ,2,4-trichlorobenzene 


40.8 


NR 


1 ,3,5-trichlorobenzene 


35 


NR 



47 



WDN:AAR:th 7/31/2000 4239-55412 NIH (DHHS) Ref. No. E-209-99/0 Express Mail No. EL#121363839US 

Date of Deposit: July 31, 2000 



1 ,2,3,4,-tetrachloroben2ene 


18.2 


R 


1 ,2,3,5-tetrachlorobenzene 


18.6 


R 


1 ,2,4,5-tetrachlorobenzene 


28.8 


R 


Pentachlorobenzene 


17.8 


R 


Hexachlorobenzene 


27.1 


R 


2-chloroaniline 


175 


NR 


3-chloroaniline 


672 


NR 


4-chloroaniline 


203 


NR 


2,3-dichloroaniline 


54 


NR 


3,4-dichloroaniline 


55 


NR 


3,5-dichloroaniline 


149 


NR 


3,4, S-trichloroaniline 


145 


NR 


2,3,4,5-trichloroaniline 


57 


NR 


Pentachloroaniline 


40 


NR 


2-chlorophenol 


6.9 


R 


3-chlorophenol 


29.2 


R 


4-chlorophenol 


11.7 


R 


2,3-dichlorophenol 


13.2 


R 


2,4-dichlorophenol 


15.3 


R 


2,5-dichlorophenol 


46.3 


NR 


2,6-dichlorophenol 


5.0 


R 


3,4-dichlorophenol 


69.9 


NR 


3,5-dichlorophenol 


17.3 


R 


2,3,4-trichlorophenol 


1.8 


R 


3,4,5-trichlorophenol 


21.7 


R 


2,3,4,5-tetrachlorophenoI 


6.5 


R 


Pentachlorophenol 


2.1 


R 
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Spectral data for these chlorobenzene compounds is obtained from the 
Integrated Spectral Data Base System for Organic Compounds web site, 
www.aist.go.jp/RIODB/SDBS, the Aldrich Library of^^C and FTNMR Spectra 
5 (Poucher and Behnke, Eds., Aldrich Chemical Company, Volumes 1-3, 1993) and the 
NIST MS database software version 1.6, Experimental ^^C NMR, EI MS, and IR 
data is collected when spectral data is not available in a database. Experimental 
spectral data is collected using standard spectroscopic protocols. 

The spectral data points are segmented and used in the same way as 3D-QSAR 

10 uses comprehensive descriptors for structural and statistical analyses (CODESSA) 
(See, long et al., /. Med. Chem., 39: 380-387, 1995 and CoUantes et al., 7. Anal 
Chem, , 68: 2038-2043, 1996, both of which are incorporated herein by reference). 
Specifically, the ^^C NMR, El-mass spectral, and IR data are placed into bins. Mass 
spectrometric data from m/z of 50 to 549 are used. Unassigned ID ^^C NMR 

15 chemical shifts are segmented into bins over a 0 to 221 ppm range. The ^^C NMR 
frequencies are shifted to bins 550 to 770, so bin 550 is the ^^C NMR spectrum for 
frequencies inside 0 to 1 ppm and bin 770 is the ^^C NMR spectrum for frequencies 
inside 220 to 221 ppm. The unassigned IR frequencies in the range from 4000 cm'^ to 
600 cm'^ are segmented into bins of 50 cm'^ width. The IR data is shifted to bins 771 

20 through 839 so that bin 771 holds the spectral data for the range 4000 cm'^ to 3950 
cm"^ and bin 839 holds the spectral data for the range 650 cm"^ to 600 cm'\ As an 
alternative, the IR data may be segmented into ranges of frequency or wavelength 
rather than the customary wavenumber frequencies. 

To save space, the ^^C NMR spectra are saved as sets of ordered pans, each 

25 consisting of the respective chemical shift frequency in ppm and the number of peaks 
within the frequency range. Likewise, the IR spectra are saved as ordered pairs of 
bin number and the number of distinct IR peaks appearing in that particular range. 
Frequency ranges with only one peak are first assigned an area of 1; doubly 
degenerate frequencies an area of 2; and so forth. This is done to provide all the 
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Spectra with similar signal-to-noise ratios and to de-emphasize line-width variations 
due to differences in experimental conditions. The bin thus defines the number of 
significant and distinct spectral features within a frequency range. Normalization of 
the NMR spectroscopic data and the IR data is accomplished by multiplying each 

5 set by a number that makes the maximum value in their respective set of bins equal to 
100, which is the maximum value for EI MS data. The number of bms used to input 
•^C NMR spectra and IR spectra may be varied to improve the SDAR. Increasing the 
number of bins and shrinking the frequency width of each bin provides separate bins 
for closely spaced spectral features that may prove important as structure descriptors 

10 for establishing a reliable SDAR. 

The segmented ^^C NMR, EI MS, and IR spectral data along with the 
biodegradability data is input into the pattern recognition program RESolve Version 
1.2 (Colorado School of Mines, Boulder, CO) and is auto-scaled and Fisher-weighted 
prior to principal component analysis (PCA). The discriminant analysis is performed 

15 based upon the canonical variate vectors. Leave-one-out (LOO) cross-validation is 
used to maximize the size of the training set. Spectral data for a test compound is 
then input and the compound judged to be readily biodegradable or not readily 
biodegradable based upon the presence of a majority of spectral features falling in 
bins that correspond to canonical variate factors from the training set that bias toward 

20 readily biodegradable or not readily biodegradable. 

Alternatively, an expert system (self-learning or not) or artificial neural 
network may be utilized to perform the pattern recognition. Where the endpoiht may 
be the result of multiple underlying mechanisms, such as multiple pathways of 
biodegradation, expert systems or artificial neural networks may be advantageously 

25 utilized to separate compounds based upon their mode of biodegradation in addition to 
their biodegradatability. 

One example of an expert system is the MuRES system that is also part of the 
RESolve 1.2 software package. The MuRES package uses spectral data that is first 
compressed by projection onto a set of principal components. This method of 
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compression is utilitarian, because it reduces the number of variables while 
maximizing the information content. MuRES uses the spectral scores that are 
calculated by projecting the spectral data onto a set of eigenvectors. However, 
MuRES may be applied directly without the compression step to data that already is 
5 overdetermined (i.e. more compounds and endpoint values than spectral data structure 
descriptors). 

The knowledge base created by MuRES is in the form of simple binary rules. 
Binary rules use binary logic, logic that can be only true or false. A complex solution 
to a problem may be decomposed into a tree, often referred to as a classification tree, 
10 which consist of simple binary rules. The largest advantage of the expert system is 
that the scores do not have to be linearly separable which is an assumption requked 
by discruninant analysis. Further details of the MuRES method are found in 
Harrington, RESolve Software Manual, Colorado School of Mines, Golden CO, 
which is incorporated herein by reference. 

15 

Example 4 - Use of NMR and UV/Visible Spectroscopic Data to Produce a 
Predictive Model for the Photosensitized Production of Singlet Molecular 

Oxygen. 

20 Phototosensitized oxidations involving singlet oxygen, a strong oxidant, are 

implicated in photodynamic inactivation of viruses and cells, in phototherapy for 
cancer, in photocarcinogenesis and m photodegradation of dyes and polymers. 
Quenching of excited singlet and triplet states of many substances by ground state 
molecular oxygen produces singlet oxygen; the lowest electronically excited singlet 

25 state of molecular oxygen. A compilation of the quantum yields for the formation of 
singlet oxygen in fluid solutions for over 700 substances is available from the Notre 
Dame Radiation Laboratory - Radiation Chemistry Data Center 
(http://www.rcdc.nd.edu). 

Compounds that are capable of photosensitizing the production of singlet 

30 oxygen are quite diverse and include aromatic hydrocarbons, aromatic ketones and 
thiones, quinones, coumarins, fluoresceins, transition metal complexes, and 
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heterocyclics. Porphyrins and pthalocyanines are particularly important classes of 
compounds that are capable of producing singlet oxygen upon illumination. 

To establish an SDAR useful for predicting the ability of molecules to 
photosensitize the production of singlet oxygen, the quantum yields for production of 

5 singlet oxygen by 20 compounds of diverse structure are obtained from the Radiation 
Chemistry Data Center website. The compounds are divided into two endpoint 
classes based upon having a high quantum yield for production of singlet oxygen (H, 
QY> 0.50) or a low quantum yield for production of singlet oxygen (L, QY< 0.50). 
Spectral data for these compounds is obtained from the resources mentioned in 

10 Examples 1 and 2 or is gleaned from other literature sources or is measured 
experimentally. The compounds, their singlet oxygen quantum yields, and the 
endpoint class for each are listed in Table 4. 



52 



WDN-AAR-th 7/31/2000 4239-55412 NIH (DHHS) Ref. No. E-209-99/0 Express Mail No. EL#121363839US 

Date of Deposit: July 31, 2000 



Table 4 



Compound 


LiuanTum Yieiu 


uriUpUlllL ^laoo 


Acetophenone 


U.zy 


1 

L 


Acridine 


0.83 


LJ 

H 


Anthracene 


0.83 


It 

hi 


9J0-Anthraquinone 


0.15 


1 

L 


Anthrone 


0.25 


L 


Azulene 


0 


L 


Benzil 


0.57 


H 


Benzophenone 


0.29 


L 


Biacetyl 


0.57 


H 


Biphenyl 


1 .0 


H 


Flourene 


0.74 


H 


1-Methylindole 


/-V OP" 

0.35 


L 


1 ,3-Diphenylisobenzofuran 


u.zo 


1 


Phenalen-1-one 


0.94 


H 


Pivalothiophenone 


1.0 


H 


Psoralen 


0.0055 


L 


Pyrene 


0.71 


H 


(E)-Stilbene 


0.08 


L 


Tetracene 


0.70 


H 


Xanthen-9-one 


0.33 


L 



Ultraviolet-visible (UV-Vis) spectral data is used along with ^^C NMR spectral 
data in the SDAR to include a measure of the importance of molecular excited states 
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to the production of singlet oxygen. The segmented spectral data is used as a set of 
descriptors in the same way that 3D-QSAR uses comprehensive descriptors for 
structural and statistical analyses. (CODESSA) (See , Tong et al., /. Med. Chem,, 
39: 380-387, 1995 and CoUantes et al., /. Anal Chem,, 68: 2038-2043, 1996, both 

5 of which are incorporated herein by reference). The ^^C NMR spectral data is 

segmented into Ippm bins over a 0 to 222ppm range. The UV-Vis data is segmented 
into 5 nm bins over the range of 190 nm to 900 nm. The ^^C NMR spectral data 
occupies bins 1 through 221 and the UV-Vis data occupies bins 222 through 321. To 
save space, the ^^C NMR spectral data is saved as sets of ordered pairs, each 

10 consisting of the bin number and the number of peaks with the frequency range 

corresponding to the bin. As described above, the ^^C NMR bins define the number 
of significant and distinct spectral features within a frequency range. The UV-Vis 
data is saved as sets of ordered pairs, each consisting of the bin number and the 
average molar absorptivity of the molecule within the wavelength range 

15 corresponding to the bin. Molar absorptivity is used rather than absorbance to correct 
for variations in concentration between the measured UV-Vis spectra of the 
compounds. 

The segmented ^^C NMR and UV-Vis spectral data along with the quantum 
yields for singlet oxygen production for 20 compounds (the traming set) are input as 

20 text files into the pattern recognition program RESolve Version 1.2 (Colorado School 
of Mines, Boulder CO). The spectroscopic data is auto-scaled and Fisher-weighted 
prior to principal component analysis (PC A). The discriminant analysis is performed 
based upon the canonical variate vectors. Leave-one-out (LOO) cross-validation is 
used to maximize the size of the training set and measure the ability of the SDAR to 

25 classify compoxmds correctly into their endpoint class. The number of bins used to 
input ^^C NMR spectral data and UV-Vis spectral data may be varied. Increasing the 
number of bins will shrink the frequency width of each bin and provide separate bins 
for closely spaced spectral features that may prove important as structure descriptors 
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for establishing a reliable SDAR. Increasing the number of compounds in the 
training set may also lead to a more reliable SDAR. 

Once the training set has been cross-validated and deemed reliable, test 
compounds are then subjected to ^^C-NMR and UV/Visible Spectroscopy, and these 

5 spectra are segmented into bins in the same manner as with the training set. Spectral 
patterns in these spectra for the test compound are then compared to the patterns for 
the training set, and endpoints associated with the spectral patterns of the training set 
are used to predict endpoints for the test compound. Where the endpoint is 
photosensitization of singlet oxygen, the presence or absence of spectral patterns in 

10 the training set which are associated with that endpoint are then detected to predict 
whether the test compound would be likely to have that characteristic. 

Example 5 - Use of an Experimentally Derived SDAR and ^^C NMR Spectral 
Data to Screen a Combinatorial Library of Compounds for Estrogen Receptor 
15 Binding Ability 

Combinatorial chemistry refers to methods of generating large numbers of 
compounds from smaller building block compounds that play an important role in 
generating lead compounds for rational drug design. The building block compounds 

20 are allowed to react to yield new compounds, either by mixed synthesis using all 
building blocks together at once or by sequential reaction of the building block 
compounds. The reactions can take place within the virtual environment of a 
computer or they can actually be carried out in a reactor system. The resulting set of 
compounds is termed a combinatorial library. 

25 A combinatorial library of compounds is produced using a combinatorial 

chemistry software package such as Afferent Structure™ (Afferent Systems, San 
Francisco, CA). With this software package it is possible to begin with a molecule 
and perform virtual synthetic steps on that molecule and any intermediate compounds 
made therefrom to yield a large number of intermediate and final product stractures. 
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For each of the structures generated using the combinatorial chemistry 
software, "C NMR spectra are predicted. The "C NMR spectra may be predicted by 
any known method. Examples of methods for predicting "C NMR spectra include 
the neural network methods described by Kvasnicka (Kvasnicka, V., J. Math. Chem., 

5 6: 63-76, 1991) and the quantum mechanical calculations of Dios et al. (Dios et al., 
Science 260:1491-1496, 1993). Software for predicting ^^C NMR spectra is also 
available from Advanced Chemistry Development, Toronto, Ontario, Canada 
(www.acdlabs.com) (ACD/CNMR Spectrum Generator). 

The predicted '^C NMR spectra are segmented into bins as the experimental 

10 '^C NMR spectra in Example 1 were segmented. The segmented, predicted spectral 
data is input into the pattern recognition program and the SDAR established for 
estrogen receptor binding based on '^C NMR only, from Example 1, is used to 
classify the compounds in the combinatorial library as either strong or medium 
estrogen receptor binders. Those compounds predicted by the SDAR to be strong 

15 estrogen receptor binders may then be tested experimentally for ability to bind to 
estrogen receptors. 

Combinatorial libraries offer a quick approach to generating large numbers of 
new compovmds, yet screening those compounds for specific biological activities is 
difficult and time consuming. The methods of the present invention provide rapid 
20 methods of screenmg the products produced by combinatorial chemistry methods so 
that time consuming assays are performed only on compounds likely to exhibit the 
desired endpoint properties. 

WhQe this embodiment utilized an SDAR generated using experimental 
spectral data, in another embodiment, calculated spectral data is used to establish the 
25 SDAR. An SDAR established with calculated spectral data may be desirable when 
screening compounds according to their calculated spectral data, as the errors in 
calculated peak position are advantageously similar in both the training set spectral 
data and test compound spectral data. 
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Example 6 - Screening of Fractions for Specific Activities using SDAR 

The SDAR methods of the present invention are useful for screening raw 

fractions of compounds derived either from natural sources or from chemical reaction 
5 mixtures (including experimental combinatorial libraries). 

A biological source of potentially new compounds, such as a sponge, is 

homogenized and partitioned to provide an aqueous fraction and an organic fraction. 

Each of these fractions is then chromatographed by any known method to yield a 

larger number of fractions that contain one or more compounds. Within some of 
10 these fractions, a compound or compounds capable of binduig to estrogen receptors 

may be present. SDAR is used to quickly screen the fractions for the presence of 

estrogen receptor binding compounds. 

An unassigned ^^C NMR spectrum is obtained for each fraction and the 

spectrum is segmented into bins of Ippm width. The fraction spectral data is then 
15 used along with the SDAR based on ^^C NMR data only, as in Example 1, to predict 

the fraction's ability to bind to estrogen receptors. Screening of fractions in this 

manner avoids the time-consuming and costly step of subjecting each fraction to a 

bioassay. 

Those fractions identified as containing compounds having the potential to 
20 bind to estrogen receptors may then be subjected to further analysis to reveal the 
number and identity of the compounds within the fraction. 

Example 7- Endpoints 

Endpoints for use with the SDAR methods encompass the full range of 
25 biological, chemical, and physical properties exhibited by molecules. The methods of 
the present invention can be used to assist in drug design, biological activity 
predictions, toxicological predictions, chemical reactivity predictions, and metabolic 
pathway predictions. An endpoint is any molecular property or activity that can be 
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measured qualitatively or quantitatively. Endpoints may be expressed in absolute or 
relative terms. 

Endpoints may be chosen to establish SD ARs that can be used to predict the 
environmental fate and toxicity of compounds. The ability of compounds to penetrate 
5 membranes, bind to enzyme active sites, react with soil, air, or water constituents, 
bind to soil constituents, hydrolyze, oxidize, and be transported in the environment 
can be used, along with spectral data for those compounds, to produce useful SDARs. 

Spectral data can be used in combination with non-specific measures of 
toxicity, mutagenticity, teratogenicity, and carcinogenicity to establish SDARs. One 

10 example of such a non-specific measure is the Ames test. DNA damage and repair 
tests, Phosphorous-32 postlabeling, mutation induction in transgenes are others. Yet 
others include transgenic mouse assays, including the p53+/- deficient model, the 
Tg. AC model, the TgHras2 model, and the XPA deficient model. LD50 and EC50 
may provide endpoints for SDAR methods as well. Alternatively, the ability of 

15 compounds to induce specific biological outcomes such as cellular changes can be 
chosen as the endpoint used to establish the SDAR. For example, relevant tissues 
may be examined for changes at the cellular level using morphological, 
histochemical, or functional criteria. As appropriate, attention may be directed to 
such changes as the dose-relationships for apoptosis, cell proliferation, liver foci of 

20 cellular alteration, or changes in intercellular communication. 

An SDAR may be established based upon any measurable response elicited in 
animals, plants, and microbes upon exposure to a series of compounds. Examples 
include SDARs based upon antiviral and antimicrobial activity. The ability of 
compounds to induce metabolic disorders such as alterations in sugar metabolism may 

25 provide a useful endpoint. Phytotoxicity and stimulation of plant growth and 
reproduction are other examples. Pesticidal activity is yet another example. 
Measures of anti-hypertensive activity, anti-pyretic activity, anti-depressant activity, 
and the like further illustrate useful endpoints that are usually related to human health. 
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Phototoxicity, both specific and non-specific, may be correlated with spectral 
features to yield an SDAR. 

Multiple endpoints may be utilized to establish multiple SDARs from a single 
set of spectral data. Compounds then may be screened based upon their spectra using 
5 multiple SDARs for any combination of desirable or undesirable activities. One 
example of a useful combination is that of maximal potential efficacy as a therapeutic 
agent with minimal potential side effects. Agrochemicals may be screened using 
multiple SDARs for species-specific toxicities and tolerances. 

An especially useful application of the methods of the present invention is to 
10 the prediction of ligand-target molecule binding. The binding of a molecule to a 
target such a protein, nucleic acid, synthetic polymer, chimeric molecule, or 
membrane constituent is often the most important step in the elicitation of a particular 
property or activity by a molecule. Binding affinities for ligand-target molecule 
interactions can be expressed in either absolute (e.g., an equilibrium constant) or in 
15 relative (e.g., relative to a reference compound, as determined for example by a 
competitive binding experiment) terms. Example 1 above is one example of how the 
relative binding affinity of a series of molecules can be utilized along with spectral 
data to establish a predictive SDAR model. SDAR models based upon relative 
binding affinities may be useful for rapidly and inexpensively screening compounds 
20 for a particular activity. They also may be useful tools for rational drug design when 
used to identify the spectral, and thus structural, features responsible for that activity. 

The metabolic pathway involved in the production or destruction of a series of 
molecules is another endpoint useful for the methods of the present invention, A 
predictive SDAR based upon pathway-structure relationships may be able to predict 
25 the biosynthetic path for newly discovered naturally occurring compounds. Similarly, 
SDAR using biodegradability as an endpoint may be useful for predicting the 
residence time of pollutants in the envuronment. 

Rates of reaction and other measures of reactivity, such as site of reaction on a 
molecular structure, including the site of electrophilic aromatic substitutions on 



59 



WDN:AAR:th 7/31/2000 4239-55412 NIH (DHHS) Ref. No. E-209-99/0 Express Mail No. EL#121363839US 

Date of Deposit: July 31, 2000 

aromatic compounds, are useful chemical endpoints for the practice of the present 
invention. 

Physical constants such as water-octanol partition coefficients, vapor 
pressures, pKa, pKb, hydrophobicities, relative acidities and basicities as well as 
5 water solubilities can be used with spectral data to provide SDARs. Such estimates 
may be especially useful for physical properties that are difficult and time-consuming 
to measure. For example, octanol- water partition coefficients are important for 
modeling the envkonmental transport of chemicals. While the octanol-water partition 
coefficient of a compound might be available, it is less likely that transient species 

10 derived from that compound during biodegradation are available m sufficient 

quantities to measure their octanol- water partition coefficients. SDAR according to 
the methods of the present invention provides an efficient way to predict the octanol- 
water partition coefficient for transient species, whose environmental transport 
characteristics need to be modeled. 

15 Other examples of endpoints useful for the methods of the present invention, 

in addition to those discussed above, may be found in Hansch and Leo, Exploring 
QSAR: Fundamentals and Applications in Chemistry and Biology, American Chemical 
Society, 1995. Further examples of endpoints useful for the methods of the present 
invention may be found in Quantitative Structure-Activity Relationships in 

20 Environmental Sciences-Vll, Chen and Schxiurmann, eds., SETAC Press, 1997. 

Example 8 - Spectral Data 

Spectroscopy refers to branch of analytical chemistry in which atomic and 
molecular structure is studied by measuring radiant energy absorbed or emitted by a 
25 substance in any of the wavelengths of the electromagnetic spectrum, in response to 
excitation by an external energy source. The types of absorption and emission 
spectroscopy are usually identified by the wavelength involved, such as gamma-ray, 
X-ray, UV, visible, infirared, microwave, and radiofrequency. Nuclear magnetic 
resonance spectroscopy (NMR) examines differences in energy states created by a 
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magnetic field. Spectral data refers to the measurements of the energy differences 
across the spectrum, and spectral patterns refer to differences in the detected energy 
differences measured across a region of the electromagnetic spectrum. Any 
instrumental method that produces data that depend upon the structural and quantum 
5 mechanical properties of a molecule may be utilized with the methods of the present 
invention. 

Spectral data as used in some embodiments includes the entire spectrum (or 
spectra) generated by the instrumental method (or methods) of spectroscopy or by 
calculation. Furthermore, the spectral data need not be assigned to particular 

10 structural features. In other embodiments the spectral data comprises only a portion 
of the spectrum or spectra available. The spectral portions utilized in the methods of 
the present invention may conveniently cover a spectral region known to typically 
arise from one or more particular structural features. For example, with respect to 
^^C NMR spectral data, spectral data can be obtained from the entire ^^C NMR 

15 spectrum (0 to 220 ppm), or at least half or a third of that spectrum, or at least a 60 
ppm, 80 ppm, 100, or 150 ppm portion of the spectrum. With respect to IR spectral 
data, for example, the spectral data can be obtained from the entire IR spectrum 
(4000 cm'^ to 500 cm"^), or at least a hundredth, fiftieth, quarter, or half of that 
spectrum, or at least a 35, 50, 100, 200, 500, or 1000 cm'^ portion of the spectrum. 

20 In particular embodiments, all or substantially all of the spectral data within a 

particular portion of the spectrum is obtained. However, the spectral features within 
the entire spectrum of within the portion of the spectrum need not be assigned to 
structural features and referenced to the corresponding spectral features arising from 
the structure features of a reference compound. For example, rather than utilizing the 

25 assigned ^^C NMR resonances of the carbon atoms in a substituted benzene ring 

referenced to the ^^C NMR resonances of the carbon atoms in unsubstituted benzene, 
the present invention could utilize the entire spectral region wherein benzene ring 
carbons typically fall, without assigning and referencing the data therein to the 
corresponding spectral data for benzene, which has a known structure. 
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Nuclear magnetic resonance (NMR) data is especially attractive for use with 
the present invention because of the large amount of structural information contained 
in a NMR spectrum. NMR instrumentation is widely available and NMR spectra are 
obtained routinely during structure elucidation. Additionally, the NMR spectra of 
5 many compounds have already been measured and are available (see Example 1 for 
some representative sources). 

^^C NMR and NMR spectral data (especially ^^C NMR spectral data) are 
very sensitive to subtle changes in substitution, conformation, chirality, and electronic 
density. Moreover, changes in ^^C NMR chemical shifts can occur at a site as many 

10 as five carbon atoms removed from the site of the variation. Solvation and proton- 
exchange effects on the electronic properties of molecules are more clearly reflected 
in NMR chemical shifts and line widths. One-dimensional ^^C NMR and NMR 
spectral data as well as two-dimensional ^^C-^H heterocorrelation data such as that 
derived from HSQC and HMQC experiments are useful. Furthermore, modern FT- 

15 NMR instruments are capable of providing NMR spectral data for as little as 1 
nanogram of a compound (an amount likely insufficient for performing standard 
bioassays such as the Ames test). 

NMR data may be segmented into bins prior to their analysis, along with 
endpoint data, in a pattern-recognition program. Suitable bin widths will vary 

20 according to the identity of the nuclei for which the spectrum is generated, and 
whether the technique is one or two-dimensional. For one-dimensional ^^C NMR 
spectral data, the bin width may be varied from the digital resolution of the 
instrument (typically about 0.1 ppm) to about 50ppm. For one-dimensional NMR 
spectral data, the bm width may be varied fi-om the mstnmiental digital resolution 

25 (typically about 0.01 ppm) to about 2 ppm. For two-dimensional ^^C-*H 

heterocorrelation data, the bin may be defined by similar corresponding widths in 
both the ^^C and dimensions. Even data of higher dimensions (e.g., three, four, 
etc.) including NMR spectral data fi^om other nuclei, such as ^^N, ^^P, ^^F, ^^O, and 
^^S may be used in establishing an SDAR. Correspondingly, bins may be defined 
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with respect to each dimension and may be of a width equal to the digital resolution 
of the data or greater. 

In anothet embodiment of the invention, ^^C NMR spectral data are predicted 
by calculation (see, for example, Dios et aL, Science 260:1491-1496, 1993 and 

5 Kvasnicka, V. , /. Math, Chem, , 6: 63-76, 1991) and used in an SDAR model that 
has been trained on true ^^C NMR spectral data. Software for predicting ^^C NMR 
spectra is also available from Advanced Chemistry Development, Toronto, Ontario, 
Canada (www.acdlabs.com) (ACD/CNMR Spectrum Generator). Predicted ^^C NMR 
spectral data may be used, for example, to aid in rational drug design by allowing 

10 proposed structures to be tested for potential activities before synthesis is attempted. 

If predicted spectral data is utilized to establish an SDAR, the spectral data 
may be segmented into bins that are desirably of a width equal to the average standard 
deviation in chemical shift predicted by the method, or greater. Likewise, test 
structure predicted spectra may be segmented in a similar fashion. 

15 Mass spectrometry can provide a measure of the size of a molecule, the size 

and identity of a molecule's structural subunits, and information regarding bond 
strengths within a molecule. Mass spectral data, especially electron impact mass 
spectral (EI MS) data, has already been obtained for many compounds and, even 
more so than NMR data, is available from convenient sources (see Example 1). EI 

20 MS data is also a standard technique used in structure elucidation. 

Other mass spectrometric techniques that are useful for providing additional 
and often complementary information include time-of-flight mass spectrometry (TOP 
MS), chemical ionization mass spectrometry (CI MS), fast-atom bombardment 
(FAB). Modern TOP MS spectrometers are capable of providing mass-spectral data 

25 from 1 ng or less of purified material (an amount that is likely insufficient for 
performing standard activity screens such as the Ames test). 

Mass spectral data may be segmented into ranges of m/z ratio (for instance, 
ranges corresponding to mass ranges from about the digital resolution of the 
instrumental method, typically about 0.1 amu, to about 50 amu) or may be segmented 
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according to integer m/z ratios, with non-integer m/z ratios being rounded to the 
nearest integer. 

Infrared (IR) spectra may also be used, for example in establishing SDAR 
models capable of discerning differences in activity seen for tautomers that are 
5 indistinguishable in NMR and mass spectrometric data. Infrared (IR) spectra are 
treated in a similar fashion to the NMR spectral data of Example 1 in that each 
spectrum may be separated into bins of a certain spectral range, for example from 
about 1 cm'^ to about 200 cm"^ for entry into a pattern-recognition program 

Ultraviolet- Visible (UV-Vis) spectral data, which is inherently reflective of the 

10 electronic energy levels of a molecule, may be used by segmenting the spectral data 
into bins of a certain spectral range, for example from about the digital resolution of 
the instrumental method to about 50 nm. UV-Vis spectral data can be used, for 
example, in predicting phototoxicity under solar illumination. Similarly, fluorescence 
and phosphorescence spectra may be handled analogously to UV-Vis spectra and 

15 utilized to establish an SDAR. Fluorescence and phosphorescence spectra reflect the 
energy redistribution within a molecule upon absorption of light and thus may provide 
important structure descriptors for predicting the light driven properties of molecules. 

In addition to the possibility of using a single type of spectral data to establish 
an SDAR, spectral data of various types may be combined to form composite sets of 

20 spectral data. Entire spectra or particular regions of spectra may be combined to 
yield spectral data sets that may be used in the methods of the present invention, 
along with endpomt data, to establish the SDAR. Spectral data may come from any 
composite of NMR, MS, IR, Fluorescence, Phosphorescence, and UV-Vis spectra, 
including composites of different species of spectra within these broad genera of 

25 spectra. Furthermore, different regions of any type of spectrum can be segmented 
into bins of different sizes so, for example, portions of a spectrum with many closely 
spaced peaks can be described by narrow spectral bins and portions of a spectrum 
without many peaks can be described by wide spectral bins. 
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In some embodiments the spectral data is not used in its raw form to establish 
an SDAR, but rather the data is subjected to pattern recognition analysis after some 
sort of pre-treatment to improve the ability of pattern recognition to extract the 
SDAR. For example normalization may be used to equalize the importance of 

5 spectral data derived from different instrumental methods when forming a composite, 
such as, a composite of MS data and NMR data wherein the maximum signals might 
be 100 and 1000 respectively. Autoscaling may be used to equalize the importance of 
inherently weak spectral data with inherently strong spectral data, for example, UV- 
Vis absorption bands within an absorption spectrum with very different extinction 

10 coefficients. Fisher-weighting may be used to emphasize the spectral data that are 
most important for predicting the endpoint data, such as spectral data found in 
compounds with a large endpoint values but absent from compounds with small 
endpoint values. 

In general, techniques for pre-treatmg data include artifact removal and/or 
15 linearization, centering, and scaling and weighting. A common form of artifact 
removal is baseline correction of a spectrum. Common linearizations include the 
conversion of spectral transmittance into spectral absorbance and the multiplicative 
scatter correction for diffuse reflectance spectra. Centering, sometimes called mean 
centering is simply the subtraction of the mean spectral signal at each frequency or 
20 m/z from each spectrum. Scaling or weighting involves multiplying all of the spectra 
by a different scaling factor for each sub-spectral region. This is done to increase or 
decrease the influence of certain spectral regions or features. A particular example of 
weighting is Fisher- weighting. Two types of scaling are typically encountered, 
variance scaling and autoscaling. Further discussion of data pretreatment may be 
25 found in Kramer, R. , Chemometric Techniques for Quantitative Analysis, Marcel 

Dekker, Inc., 1998. Additional methods for the pre-treatment of data prior to pattern 
recognition are known in the art and are mcluded in the present invention. 

Since spectral data is utilized in some embodunents in the same manner as 3D- 
QSAR utilizes comprehensive descriptors for strucmral and statistical analyses 
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(CODESSA) (see, Tong et al., /. Med, Chem., 39: 380-387, 1995 and CoUantes et 
aL, 7. Anal Chem, , 68: 2038-2043, 1996, both of which are incorporated herein by 
reference) it is possible to combine the spectrally derived structure descriptors of the 
present invention v^ith structure descriptors based on bulk measurements (such as the 

5 octanol/water partition coefficient) and still retain the ability to produce predictive 
models without the need for structural knowledge beforehand. Hybrid methods also 
fall within the scope of the invention wherein the spectrally derived descriptors, 
which do not require structural knowledge beforehand, may be combined with other 
structure descriptors that require structural knowledge to produce a larger set of 

10 descriptors for use in a predictive model. 

Example 9 - Pattern-Recognition Programs 

Pattern-recognition programs useful for practicing the present invention are of 
two major types; statistical and artificial intelligence. 

15 Statistical methods include Principal Component Analysis (PCA) and 

variations of PCA such as linear regression analysis, cluster analysis, canonical 
variates, and discriminant analysis, soft independent models of class analogy 
(SIMCA), expert systems, and auto spin (see, for example, Harrington, RESolve 
Software Manual, Colorado School of Mines, 1988, incorporated by reference). 

20 Other examples of statistical analysis software available for principal-component- 
based methods include SPSS (SPSS Inc., Chicago, IL), JMP (SAS Inc., Cary NC), 
Stata (Stata Inc., College Station, TX) and Cluster (available to run from 
entropy: ~dblank/public_htnil/cluster). 

Artificial intelligence methods include neural networks and fuzzy logic. 

25 Neural networks may be one-layer or multilayer in architecture (See, for example, 
Zupan and Gasteiger, Neural Networks for Chemists, VCH, 1993, incorporated herein 
by reference). Examples of one-layer networks include Hopfield networks. Adaptive 
Bidirectional Associative Memory (ABAM), and Kohonen Networks. Examples of 
Multilayer Networks include those that learn by counter-propagation and back- 
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propagation of error. Artificial neural network software is available from, among 
other sources, Neurodimension, Inc., Gainsville, FL (Neurosolutions) and The 
Mathworks, Inc., Natick, MA (MATLAB Neural Network Toolbox). 

Spectral patterns can be analyzed using other approaches. Analog spectral 
5 peak patterns may be digitized, and image analysis may used to search for similarities 
or differences between the spectral patterns of training sets and test compounds 

Example 10 - Computer Environments for SDAR Methods. 

The SDAR methods of the present invention may be implemented using a 

10 single computer or utilizing a distributed computing environment. 

Exemplary Distributed Computing Environment 
Fig. 13 illustrates a distributed computing environment in which the software 
elements used to implement the SDAR methods of the present invention may reside. 
The distributed computing environment 100 includes two computer systems 102, 104 

15 connected by a connection medium 106. The computer systems 102, 104 can be any 
of several types of computer system configurations, including personal computers, 
multiprocessor systems, and the like. In terms of logical relation with other computer 
systems, a computer system can be a client, a server, a router, a peer device, or other 
common network node. Moreover, although Fig. 13 illustrates two computer systems 

20 102, 104, the present invention is equally applicable to an arbitrary, larger number of 
computer systems connected by the connection medium 106. Additional computer 
systems 102 or 104 may be connected by an arbitrary number of connection mediums 
106. The connection medium 106 can comprise any local area network (LAN), wide 
area network (WAN), or other computer network, including but not limited to 

25 Ethernets, enterprise-wide computer networks, intranets and the Internet. 

Portions of the SDAR software can be implemented in a smgle computer 
system 102 or 104, with the application later distributed to other computer systems 
102, 104 in the distributed computing environment 100. Portions of the SDAR 
software may also be practiced in a distributed computing environment 100 where 
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tasks are performed by a single computer system 102 or 104 acting as a remote 
processing device that is accessed through a communications network, with the 
distributed application later distributed to other computer systems in the distributed 
computing environment 100. In a networked environment, program modules 
5 comprising the SDAR software can be located on more than one computer system 102 
or 104. Communication between the computer systems in the distributed computing 
network may advantageously include encryption of the communicated data. 

Exemplary Computer System 

10 Fig. 10 illustrates an example of a computer system 120 that can serve as an 

operating environment for the SDAR software. With reference to Fig. 14 an 
exemplary computer system for implementing the invention includes a computer 120 
(such as a personal computer, laptop, palmtop, set-top, server, mainframe, and other 
varieties of computer), including a processing unit 121, a system memory 122, and a 

15 system bus 123 that couples various system components including the system memory 
to the processing unit 121. The processing unit can be any of various conmiercially 
available processors, mcluding Intel x86, Pentium and compatible microprocessors 
from Intel and others, including Cyrix, AMD and Nexgen; Alpha from Digital; MIPS 
from MIPS Technology, NEC, IDT, Siemens, and others; and the PowerPC from 

20 IBM and Motorola. Dual microprocessors and other multi-processor architectures 
also can be used as the processing unit 121. 

The system bus can be any of several types of bus structure including a 
memory bus or memory controller, a peripheral bus, and a local bus using any of a 
variety of conventional bus architectures such as PCI, VESA, AGP, MicroChannel, 

25 ISA and EISA, to name a few. The system memory includes read only memory 
(ROM) 124 and random access memory (RAM) 125. A basic input/output system 
(BIOS), containing the basic routines that help to transfer information between 
elements within the computer 120, such as during start-up, is stored in ROM 124. 
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The computer 120 further includes a hard disk drive 127, a magnetic disk 
drive 128, e.g., to read from or write to a removable disk 129, and an optical disk 
drive 130, e.g., for reading a CD-ROM disk 131 or to read from or write to other 
optical media. The hard disk drive 127, magnetic disk drive 128, and optical disk 
5 drive 130 are connected to the system bus 123 by a hard disk drive interface 132, a 
magnetic disk drive interface 133, and an optical drive interface 134, respectively. 
The drives and their associated computer-readable media provide nonvolatile storage 
of data, data structures, computer-executable instructions, etc. for the computer 120. 
Although the description of computer-readable media above refers to a hard disk, a 
10 removable magnetic disk and a CD, it should be appreciated by those skilled in the art 
that other types of media which are readable by a computer, such as magnetic 
cassettes, flash memory cards, digital video disks, Bernoulli cartridges, and the like, 
can also be used in the exemplary operating environment. 

A number of the SDAR program modules can be stored in the drives and 
15 RAM 125, including an operating system 135, one or more application programs 136, 
other program modules 137, and program data 138. 

A user can enter commands and information into the computer 120 through a 
keyboard 140 and pomting device, such as a mouse 142. Other input devices (not 
shown) can include a microphone, satellite dish, scanner, or the like. These and 
20 other input devices are often connected to the processmg unit 121 through a serial 
port interface 146 that is coupled to the system bus, but can be connected by other 
interfaces, such as a parallel port, game port or a universal serial bus (USB). A 
monitor 147 or other type of display device is also connected to the system bus 123 
via an interface, such as a video adapter 148. In addition to the monitor, computers 
25 typically include other peripheral output devices (not shown), such as printers. 

The computer 120 can operate in a networked environment using logical 
connections to one or more other computer systems, such as computer 102. The 
other computer systems can be servers, routers, peer devices or other common 
network nodes, and typically include many or all of the elements described relative to 
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the computer 120, although only a memory storage device 149 has been illustrated in 
Fig. 14. The logical connections depicted in Fig. 10 include a local area network 
(LAN) 151 and a wide area network (WAN) 152. Such networking environments are 
common in offices, enterprise-wide computer networks, intranets and the Internet. 

5 When used in a LAN networking environment, the computer 120 is connected 

to the local network 151 through a network interface or adapter 153. When used in a 
WAN networking environment, the computer 120 typically includes a modem 154 or 
other means for establishing communications (e.g., via the LAN 151 and a gateway 
or proxy server 155) over the wide area network 152, such as the Internet. The 

10 modem 154, which can be internal or external, is connected to the system bus 123 via 
the serial port interface 146. In a networked environment, program modules depicted 
relative to the computer 120, or portions thereof, can be stored in the remote memory 
storage device. It will be appreciated that the network connections shown are 
exemplary and other means of establishing a communications link between the 

15 computer systems (including an Ethernet card, ISDN termmal adapter, ADSL 

modem, lOBaseT adapter, 100BaseT adapter, ATM adapter, or the like) can be used. 

In accordance with the practices of persons skilled in the art of computer 
programming, a particular embodiment of the SDAR method is described in Fig. 4 
with reference to acts and symbolic representations of operations that may be 

20 performed by the computer 120. Such acts and operations are sometimes referred to 
as being computer-executed. It will be appreciated that the acts and symbolically 
represented operations include the manipulation by the processing unit 121 of 
electrical signals representing data bits which causes a resulting transformation or 
reduction of the electrical signal representation, and the maintenance of data bits at 

25 memory locations in the memory system (including the system memory 122, hard 
drive 127, floppy disks 129, and CD-ROM 131) to thereby reconfigure or otherwise 
alter the computer system's operation, as well as other processing of signals. The 
memory locations where data bits are maintained are physical locations that have 
particular electrical, magnetic, or optical properties corresponding to the data bits. 
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In view of the many possible embodiments to which the principles of our 
invention may be applied, it should be recognized that the illustrated embodiments are 
only specific examples of the invention and should not be taken as a limitation on the 
scope of the invention. Rather, the scope of the invention is defined by the following 
claims. We therefore claim as our invention all that comes within the scope and spirit 
of these claims. 
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We Claim: 

A method of predicting a biological activity of a molecule, comprising: 
obtaining spectral data for a test compound; 
5 comparing the spectral data of the test compound to a pattern of spectral data 

associated with a biological activity, derived not exclusively from the assigned 
spectral data of a training set of compounds having a known biological activity; 

detecting similarities between the pattern of spectral data associated with a 
biological activity of the training set and a pattern of spectral data for the test 
10 compound to determine whether the test compound is predicted to share the biological 
activity. 

2. The method of claim 1 wherein the spectral data are obtained without 
first correlating the spectral data with corresponding structural features. 

3 The method of claim 1 wherein the pattern of spectral data associated 
15 with a biological activity is derived without first correlating the spectral data with 
corresponding structural features, 

4. The method of claim 1, wherein the pattern of spectral data of the 
training set is a pattern obtained by separating the spectral data of the training set of 
compounds into sub-spectral units. 
20 5. The method of claim 4, wherein the pattern of spectral data of the test 

compound is obtained by separating the spectral data of the test compound into 
substantially the same sub-spectral units into which the spectral data of the trainmg set 
is separated. 

6. The method of claim 1, wherein the spectral data is one type of 
25 spectral data. 

7. The method of claim 6, wherein the spectral data comprises one of 
nuclear magnetic resonance, mass spectral, infi^ared, ultraviolet-visible, fluorescence, 
or phosphorescence data. 
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8. The method of claim 1, wherein the spectral data is a composite of 
different types of spectral data. 

9. The method of claim 8, wherein the different types of spectral data 
comprise two or more of the group consisting of nuclear magnetic spectroscopy 

5 (NMR), mass spectroscopy (MS), infrared (IR) spectroscopy, and ultraviolet-visible 
(UV-Vis) spectroscopy. 

10. The method of claim 1, wherein the spectral pattern of the test 
compound and the spectral pattern of the training set are segmented into sub-spectral 
units, and the spectral data of the training set is scaled to normalize the importance of 

10 different signals within the spectral data of the training set prior to deriving a pattern 
associated with a biological activity. 

1 1 . The method of claim 10, wherein the scaling is auto-scaling. 

12. The method of claim 10, wherein the spectral data of the training set is 
weighted to emphasize signals that are important for determining the endpoint class of 

15 compounds in the training set before deriving a pattern associated with a biological 
activity. 

13. The method of claim 12 wherein the weighting is Fisher- weighting. 

14. The method of claim 1, wherein detecting similarities between the 
pattern of spectral data associated with a biological activity of the training set and the 

20 pattern of spectral data for the test compound comprises performing computer 
implemented pattern recognition. 

15. The method of claim 1, wherein detecting similarities between the 
pattern of spectral data associated with a biological activity of the training set and the 
pattern of spectral data for the test compound comprises detecting relative intensities 

25 of signals associated with one or more of the sub-units of the spectrum of the training 
set, and detecting relative intensities of signals associated with the same one or more 
sub-units of a spectrum of the test compound. 
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10 



15 



20 



16. The method of claim 15, wherein the relative intensities are canonical 
variate factors of the spectral data associated with a biological activity of the training 
set and the spectral signals of the test compound, 

17, The method of claim 1, wherein the method is computer implemented. 



a test compound, comprising: 

receiving as input spectral data for a test compound; 

receiving as input spectral data of a training set of compounds having a known 
biological activity; and 

comparing the pattern of spectral data of the training set associated with the 
biological activity to the spectral pattern of the test compound to determine whether 
the spectral pattern of the test compound matches the spectral pattern associated with 
the biological activity of the training set. 

19. The computer implemented system of claim 18, wherein comparing the 
spectral patterns comprises comparing the spectral patterns with computer 
implemented pattern recognition programs. 

20. The computer implemented system of claim 19, wherein the spectral 
data for the test compound and the spectral data for the training set are divided into 
substantially identical spectral bias, so that a signal within individual spectral bins is 
compared between spectral patterns of the training set associated with the biological 
activity and the test compound. 

21. The computer implemented system of claim 18, wherein the spectral 
patterns are obtained by inputting spectral data selected from the group consisting of 
nuclear magnetic resonance data, mass spectral data, infrared data, untraviolet-visible 
data, fluorescence data, phosphorescence data, and composites of two or more such 
spectral data. 

22. The computer implemented system of claim 21, wherein the spectral 
data for the training set are converted into canonical variates associated with the 




A computer implemented system for predicting a biological activity of 
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biological activity of the training set, and the spectral data for the test compound are 
compared to the canonical variates of the training set spectral data. 

23 . The computer implemented system of claim 22, wherein the biological 
activity is binding affinity to a hormone receptor, and the canonical variates for the 
training set include peaks in bms that are associated with hormone receptor binding of 

a pre-selected affinity. 

24. The computer implemented system of claim 23 , wherein the spectral 
data comprise nuclear magnetic resonance data and mass spectral data. 

25. A computer readable medium having stored thereon instructions for 
performing the actions of claim 1. 

26. A computer readable medium having stored thereon instructions for 
performing the actions of claim 18. 



training compounds; 

inputting the segmented spectral data and endpoint data into a pattern- 
recognition program; 

training the pattern-recognition program with the segmented spectral data and 
endpomt data to establish a relationship between the spectral sub-units of the 
segmented spectral data and the endpoint; 

providing segmented spectral data for a test compound that is segmented into 
substantially the same spectral sub-units as the spectral data of the trammg 
compounds; and 

comparmg the relationship between the spectral subunits of the segmented 
spectral data and the endpoint to the spectral subunits of the test compound's 
segmented spectral data to predict the endpoint of the test compound. 

28 . The method of claim 27, wherein knowledge of the structures of the 
training compounds and the test compound are not necessarily known beforehand. 
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29. The method of clahn 27 wherein the segmented spectral data of the 
training set is autoscaled and Fisher-weighted before training the pattern recognition 
program. 

30. The method of claim 27 wherein the spectral data is chosen from the 
5 group consisting of nuclear magnetic resonance data, mass spectral data, infrared 

data, UV-Vis data, fluorescence data, phosphorescence data, and composites thereof. 

31. The method of claim 30 wherein the spectral data is chosen from the 
group consisting of ^^C NMR data, EI MS data, and composites thereof. 

32. The method of claim 27, wherein the endpoint is a ligand-target 
10 molecule-binding affinity. 

33. The method of claim 32, wherein the ligand-target molecule binding 
affinity is an estrogen-receptor binding affinity. 

34. The method of claim 27, wherein the endpoint is selected from the 
group consisting of: 

1 5 a measure of biodegradability ; 

a measure of toxicity; 

participation in a metabolic pathway; 

a partition coefficient; 

a reaction rate; 
20 a quantum yield; 

a measure of phototoxicity; 

an equilibrium constant; and 

a site of reaction on a molecular structure. 

35. The method of claim 34, wherein the partition coefficient is the 
25 octanol-water partition coefficient. 

36. The method of claim 27, wherein non-spectral structure descriptors that 
do not necessarily depend upon structural knowledge beforehand are provided for the 
test compound and the training compounds; used to establish a relationship between 
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the segmented spectral data, the non-spectral structure descriptors, and the endpoint 
for the training compounds; and used to predict the endpoint of the test compound. 

37. The method of claim 36, wherein the non-spectral structure descriptors 
are chosen from the group consisting of partition coefficients, solubilities, relative 

5 acidities, relative basicities, pKa, pKb, reaction rates, and equilibrium constants. 

38. The method of claim 37, wherein the partition coefficient is the 
octanol-water partition coefficient. 

39. The method of claim 27, wherein non-spectral structure descriptors that 
do depend upon structural knowledge beforehand are provided for the test compound 

10 and the training compounds; used to establish a relationship between the segmented 
spectral data, the non-spectral structure descriptors, and the endpoint for the training 
compounds; and used to predict the endpoint of the test compound. 

40. The method of claim 39, wherein the non-spectral descriptors are 
calculated using a quantum mechanical or electrostatic potential method. 

15 A method for using spectral data as a set of structure descriptors for a 

compound that does not necessarily require knowledge of the compound's structure 
beforehand, comprising: 

providing spectral data; and 

segmenting the spectral data into bins. 



biological, chemical, or physical property, comprising: 
providing spectral data; 
segmenting the spectral data into bins. 

detecting patterns in the bins of the spectral data that are associated with the 



43. The method of claim 42, further comprising detecting corresponding 
patterns in spectral data of test compounds to select test compounds having the 
property. 



20 




A method for establishing a relationship between spectral data and a 



25 property. 
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44. The method of claim 43, wherein the test compounds are mixtures of 
compounds. 

45. The method of claim 42 including one or more of: auto-scaling the 
segmented data; and weighting the segmented spectral data. 

46. The method of claim 45 wherein weighting of the segmented data 
comprises Fisher-weighting of the segmented spectral data. 

47. The method of claim 1, wherein the biological activity of the test 
compound is predicted without reference to a chemical structure of the test 
compound. 



comprising; 

providing endpoint data for a plurality of compounds; 
providing spectral data for a plurality of compounds; 
segmenting the spectral data for the plurality of compounds into bins; 
autoscaling the numerical data obtained from the spectral features within each 
of the bins; 

Fisher- weighting the data within each of the bins; and 
correlating information in the bins and the endpoint using a means for pattern 
recognition. 

49. The method of claim 47 wherein the spectral data is selected from the 
group consisting of NMR data, MS data, IR data, UV-Vis data, fluorescence data, 
phosphorescence data, and composites thereof. 

50. The method of claim 48 wherein two or more types of spectral data are 
normalized to each other in a composite. 

p\. A method for determining the structural features of a plurality of 



compounds that contribute to determining a particular endpoint property exhibited by 
the compounds, comprising: 

providing segmented spectral data for the plurality of compounds; 

providing endpoint data for the plurality of compounds; 




A method for establishing a spectral data activity relationship, 
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establishing a spectral data-activity relationship by identifying the segmented 
spectral features that bias toward increased endpoint values and the segmented 
spectral features that bias toward decreased endpoint values; and 

identifying the structural feature leading to the segmented spectral features that 
5 bias toward increased or decreased endpoint values for the plurality of compounds. 

52. The method of claim 51 wherein the segmented spectral data is selected 
from the group consisting of NMR data, MS data, IR data, UV-Vis data, fluorescence 
data, phosphorescence data, and composites thereof. 

53. The method of claim 52 wherein two or more types of segmented 
10 spectral data are normalized to each other in a composite. 

54. The method of claim 52 wherein the segmented spectral data is a 
composite of ^^C NMR and EI MS data. 
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Abstract 

Methods are disclosed for establishing a quantitative relationship between 
spectral properties of molecules and a biological, chemical, or physical endpoint of 

5 the molecules. Spectral data including data from nuclear magnetic resonance, mass 
spectrometric, infrared, and ultraviolet-visible techniques are used along with 
endpoint data to train a pattern-recognition program. The training yields a spectral 
data-activity relationship that may be used to predict the endpoint value of a molecule 
from its spectral data alone. Methods for rapidly screening isolated compounds or 

10 mixtures of compounds based upon their spectral data are included. 
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FIG. 3 - a hypothetical set of spectrally derived molecular structure 

descriptors for bisphenol A. 
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Spectral Data 
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Segmentation into Bins 22 



Segmented Spectral Data 
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Auto-scaling 26 



Auto-scaled, segmented spectral data that equalizes the importance of inherently 
weak spectral signals and inherently strong spectral signals. 
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Fisher-Weighting 30 



Fisher-weighted, auto-scaled, segmented spectral data that emphasizes those spectral 
characteristics that are most important for discerning the endpoint classes of the 
compounds in the training set. 



32 



Pattern Recognition by, for example, 
Principal Component Analysis (PCA) 



The spectral data-activity relationship (SDAR), a quantitative relationship 
between segmented spectral data and endpoint data. The SDAR may be 
visualized with a canonical variate function, that illustrates how spectral data in a 
particular bin shows bias toward a particular endpoint property and a 
discruninant function that illustrates how well the SDAR is able to discern the 
endpoint group of a compound from hs segmented spectral data. 
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International appHcation in the manner provided by the first paragraph of Title 35, United States Code, § 112, 1 
acknowledge the duty to disclose material information as defined in Title 37, Code of Federal Regulations, § 1.56(a) 
which occurred between the filing date of the prior application and the national or PCT International filing date of 
this application: 



09/496,314 


February 1,2000 


Pending 


(Application No.) 


(Filing Date) 


(Status: patented, pending, abandoned) 



The undersigned hereby authorizes the U.S. attorney or agent named herein to accept and follow 

instructions from as to any action to be taken in the Patent and Trademark Office regarding this application 

without direct communication between the U.S. attorney or agent and the undersigned. In the event of a change in 
the persons from whom instructions may be taken, the U.S. attorney or agent named herein will be so notified by the 
undersigned. 

I hereby appoint the following attomey(s) and/or agent(s) to prosecute this application, to file a 
corresponding international application, and to transact all business in the Patent and Trademark Office connected 
therewith: 



Name 


Reg. No. 


Name 


Reg. No. 


James C. Haight 


25,588 


Jack Spiegel 


34,477 


Robert Benson 


33,612 


David R. Sadowski 


32,808 


Susan S. Rucker 


35,762 


Steven Ferguson 


38,448 


John Peter Kim 


38,514 


Joseph K. Hemby 


42,652 


Stephen L. Finley 


36,357 


Marlene Shinn 


46,005 


Richard R. Rodriguez 


45,980 






.ssociate Power of Attorney to the following: 






Name 


Reg. No. 


Name 


Reg. No. 


Richard J. Policy 


28,107 


Lisa M. Caldwell 


41,653 


William D. Noonan 


30,878 


Michael D. Jones 


41,879 


Mark L. Becker 


31,325 


Tanya M. Harding 


42,630 


Donald L. Stephens Jr. 


34,022 


Gregory L. Maurer 


43,781 


Stacey C. Slater 


36,011 


Paula A, DeGrandis 


43,581 


Michelle L. Johnson 


36,352 


Samuel E.George 


44,119 


Susan Alpert Siegel 


43,121 







all of the law firm of Klarquist Sparkman Campbell Leigh & Whinston, LLP. 

Address all telephone calls to William D. Noonan, M.D., telephone number 503/226-7391 and facsimile 
number 503/228-9446. 

Address all correspondence to: 

KLARQUIST SPARKMAN CAMPBELL 

LEIGH & WHINSTON, LLP 

One World Trade Center, Suite 1600 

121 SW Sahnon Street 

Portland, OR 97204-2988 
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I hereby declare that all statements made herein of my own knowledge are true and that all statements 
made on information and belief are believed to be tme; and further that these statements were made with the 
knowledge that willful false statements and the like so made are punishable by fine or imprisonment, or both, under 
Section 1001 of Title 18 of the United States Code and that such willful false statements may jeopardize the validity 
of the application or any patent issued thereon. 



Full Name of First Inventor: ^^,_J)wight W. Miller 

Inventor's Signature ^-^^-^D 

X Date 
Residence: Pine Bluff, Arkansas 

Citizensliip: United States 

Post Office Address: 2901 Dan Road 

Pine Bluff, AR 71602 

Full Name of Second Joint Inventor: Richard Beger 

Inventor's Signature /^..gy^l^^tf^ ^ /2,'//oc 



Residence: White Hail, Arkansas 
Citizenship: 

Post Office Address: 405 Reynolds #10 

White Hall, AR 71602 



Residence: Little Rock, Artcansas 
Citizenship: 

Post Office Address: 6 Pleasant Coiirt 

Little Rock, AR 72211 



Full Name of Third Joint Inventor: Jackson O. Lay, Jr. 

Inventor's Signature L M^ (L ^ , 7^6/^^ 
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Full Name of Fourth Joint Inventor: Jon G. Wilkes 

Inventor's Signature C k/| lA/jL.^^ 

Residence: Little Rock, Arkansas 
Citizenship; 

Post Office Address: 621 Cumberland #1 

Little Rock, AR 72202 

Full Name of Fifth Joint Inventor: James P. Freeman 

Inventor's Signature '^^^^''1^^^ ^ ^ ^ ^ ^ 

Residence: White Hall, Arkansas 
Citizenship: 

Post Office Address: 1 100 Rock 

White Hall, AR 71602 



Date 
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